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SYNOPSIS 

"'Experimental and theoretical investigations are made 
to study perfoimance of a gaseous rocket motor. 

Experimental work is carried out at different mass 
flow rates, at various oxidizer fuel ratios,- using nozzles of 
varying expansion ratios. The flame-in limits are also es- 
tablished. 

Theoretical investigations include (a) Determination 
of equilibrium temperature and composition of chemical reaction, 
(b) performance paraneter calculations by frozen flow technique 
and (c) perfoimance parameter calculations by equilibrium flow 
technique. \f 

Computer programmes are developed for use on IBM 7044/ 
1401 system. First programme takes care of part (a). It allows 
simultaneous high speed calculations of equilibrium composition 
and temperature for Hg - Og - Kg system, for combustion at 
constant pressure and adiabatic conditions. Second and third 
programme dealnwith part (b) and (c). They incorporate basic 
assumptions of analysis. Satisfactory results are obtained 
within maximum error limit of + 3 %• 



SECTION I 


EXPERIMENT 11 INVESTIGATIONS OE STATIC 
PERFORMANCE OE A CHEMICAL ROCKET MOTOR 



CH/PTER I 


HIT RODCJCT ION 


1.1 General 

1.2 Work Related to Performance Analysis 


1.3 Present Work 
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1.1 General 

A rocket engine is capable of operating independently 
of its environmental medium because it carries its own prope- 
llants, both fuel and oxidiser. According to the type of the 
energy source used, rocket engine can be classified as chemical 
rocket engines, nuclear rocket engines, electrical rocket en- 
gines and solar rocket engines. 

In chemical rocket engines, the energy from a high 
pressure combustion reaction of fuel and oxidising chemicals 
(propellants) adds to the heating of reaction product gases to 
a very high temperature. The high temperature gases thus ob- 
tained are subsequently expanded In a nozzle which permits con- 
version of their heat energy into kinetic energy. Physical 
state of the propellants used, once again classifies the che- 
mical rocket engines, viz. solid propellant rocket engines, 
gaseous propellant rocket engines and hybrid rocket engines. 
Controlability of liquid propellant rocket engines has made 
them more versatile. Though solid propellant rockets are used 
sometimes in missiles and RATO applications. Gaseous propell- 
ant rockets have found limited applications as attitude con- 
trol systems in space ■vehicles. Chemical additives for 
electron supression in rocket flames have been studied by 
Farber & Srivastava (l, 2). 
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1*2 Work Related to Performance .Analysis 

A great deal of work was proposed and had been suc- 
cesively carried out during and after 2nd world war. 

Sutton (3) and Barrare (4) have discussed at length 
the subject of rocket testing and many other experimental 
techniques in their text books, Sutton (5) and Young quist (6) 
also tried to represent the thermo chemistry of rocket prope- 
llants and testing of liquid propellant rocket motor as a sub- 
ject. They carried out some experiments with gasoline - 
liquid oxygen rockets and nitric acid - aniline rockets and 
after determining theoretical performance of the same, a fruit- 
ful comparision was established. 

Tenebaum (7) and Anderten (8) discussed at length 
about testing of the giant rocket engines. G-eothert (9) end 
EMI (10) also discussed the testing of rocket motor under simu- 
lated conditions. 

Mott ram (ll) made a systematic survey of measurements 
of performance of liquid rocket motors. He aLso reviewed the 
special requirements of all necessary instruments used for test- 
ing of rocket motors. Bierlin (12) summarised the methods of 
thrust measurements, both direct and indirect. Venn (l3), 

Jones (l4) and Atman (l5) worked on instrumentations for rocket 
piotor testing e.g. on various oscilloscopes, pressure recor- 
ders etc. Bowers (16) end Shafer (17) did remarkable work in 
flow measurement and calibration of flow measuring instruments* 
Exhaust gas analysis was done by Bear (18). 
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1*3 Present Work 

The present work includes the study of static perfor- 
mance of H 2 - 0 2 gaseous rocket motor unit. Poliowing major 
experiments are performed j 

(1) Injector flow analysis 

(2) Plane-in limits establishment 

(3) Thrust meter analysis 

(4) Determination of performance parameters at diffe- 
rent mass flow rates at various oxidizer fuel 
ratios using nozzles of varying expansion ratios. 



CHAPTER II 


EXPERIMENTAL 


2.1 DESCRIPTION OP APPARATUS 

2.1.1 Thrust Chamber Assembly 

2.1.2 Propellant Peed System 

2.1.3 Valves and Controls 

2.1.4 Other Components and Mounts 

2.2 FABRICATION 

2.3 EXPERIMENTAL PROCEDURE 

2.3.1 Injector Plow Analysis 

2.3.2 Plame-in Limils Establishment 

2.3.3 Thrust Meter Analysis 

2.3.4 Determination of Performance Parameters 



2.1 Description of Apparatus 
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All the experiments are performed on a IABROC VI 
assembly provided by Astrosystem International Inc. U.S.A. 

Front and rear view of the unit are shown in Fig. (l) and 
Fig. (2). This rocket motor operates on hydrogen oxygen 
(commercial gas) propellants. The assembly can be examined in 
detail with following subdivisions. 

2.1.1 Thrust Chamber Assembly 

The entire unit (see Fig. 3) is made up of high tem- 
perature nickle alloy. It consists of injector (Fig. 4), com- 
bustion chamber (Fig. 5) and exhaust nozzle (Fig. 6). 

Injector is of direct impingment type. Two jets of 
hydrogen gas are made to impinge central horizontal jet of oxy- 
gen gas at an angle of approximately 20°. Prior to as well as 
after completion of every firing, purging of combustion chamber 
is done with nitrogen gas. The combustion chamber is cylindri- 
cal, with provision for oxygen inlet, nitrogen inlet, hydrogen 
inlet, thrust knob, water inlet and outlet connections and nozzle 
mounting arrangements. Nozzles of various expansion ratios can 
be attached to the combustion chanbcr with metal *0' ring seal. 
The entire assembly is water cooled. 

2.1.2 Propellant Feed System 

Commercial gas cylinders for hydrogen, oxygen, and ni- 
trogen serves the purpose of propellant stor^e tanks. All these 
gases are made available to the combustion chamber via regulator. 









RG.4_ IMPINGING STREAM TYPE , TRIPLE 
HOLE INJECTOR 
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FIG.6_ ROCKET EXHAUST NOZZLES 












This regulates both pressure and mass flow rate of the gases. 
After regulators there follows the solonoid operated valves. 
For details see Fig. ( 7 ). 

2.1.3 Valves and Controls 

Control circuit diagram is shown in Fig. (8). 
Solenoid operated valves and other check valves are 
so connected in control circuit that, when water level is op- 
timum which puts water switch on, both oxidiser and fuel 
solenoid valves are on open position, with fire switch on. 
When shut down switch is on only nitrogen gas will purge th- 
rough gas circuits. Check switches and indicator bulbs are 
also provided on the console to check the master power, water 
level sufficiency, spark check and running time meter opera- 
tion. Water switch is adjusted to minimum level of 7 lb/min. 
flow while running time meter is set at minimum of 75 Psig. 
Belays and timers are also used in the circuit which allow 
for a purging time of 17 seconds for nitrogen gas. 

2.1.4 Other Components and Mounts 

This includes cooling water circulation pump, roto- 
meter for flow measurement, thermometers for water inlet and 
outlet temperature to the thrust chanber, chamber pressure 
measurement gauge, propellant storage tank pressure gauge, 


thrust meter, 















Nomenclature For Figure ( 7 ) 


1 Pressure regulator nitrogen purge 

2 Pressure regulator oxygen supply 

3 Pressure regulator hydrogen supply 

4 High pressure line oxygen supply 

5 High pressure line hydrogen supply 

6 High pressure line nitrogen supply 

7 Pressure gauge oxygen tank 

8 Pressure g auge hydrogen tank 

9 Pressure gauge nitrogen tank 

10 Pressure gauge oxygen supply 

1 1 Pressure g auge hydrogen supply 

12 Pressure gauge nitrogen supply 

13 Solenoid valve oxygen propellant valve 

14 Solenoid valve hydrogen propellant valve 

15 Solenoid valve nitrogen purge 

16 Check valves nitrogen purge, oxygen system 

17 Cheek valves nitrogen purge, hydrogen system 

18 Check valve oxygen system 

19 Check valve hydrogen system 

20 Pressure gauge chamber pressure 

21 Thermometer cooling water outlet temperature 

22 Thermometer cooling water inlet temperature 

23 Flow meter, cooling water 

24 Orifice, cooling water flow orifice 

25 Pressure switch, cooling water 

26 Pressure switch, chamber pressure 
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Nomenclature for Figure ( 8 ) 


F^ Fuse 

Pilot light, Power (vtoite) 

I >2 Pilot light, Purge (green) 

Pilot light, Coolant (blue) 

£4 Pilot light, firing (red) 

Ij^ Pilot light, spark check (amber) 

Mg Purge off momentary contact button 

Fire momentary contact button 
Shutdown momentary contact button 
PS.j Coolant pressure switch 

PBg Chamber pressure switch 

Purge relay UPST 

Eg Fire relsy TPIffi 

SW^ Master switch DPST 

SWg Emergency shut down switch EP32T 

SW^ Circuit spark check switch 

Transformer 

BTM Sunning time meter 

Spark Coil 
TE 


SPST 

SPST 

SPST 


Purge timer 
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2.2 Fabrication 

This includes fabrication and calibration setup of 
thrust meter. The thrust meter was fabricated to measure the 
thrust of the rocket motor. See Fig. (9). 

Bounded wire type of resistance strain gauge was 
used and it was mounted on thin steel plate of * I’ shape. The 
thrust knob from the rocket motor was made to touch the other 
side of * I' .shaped plate with the help of initial load see 
Fig. (10). 

Output cf the strain gauge was taken to strain gauge 
indicator. Proper bridge connections are essential. Calibra- 
tion of this thrust meter was done before each set of firing, 
with the help of known weights. 

2.3 Experimental Procedure 

Prior to actual test run calibration of injector for 
gas flow is done together with thrust meter calibration with 
known weights. Experiments of establishment of flame-in limits 
and performance parameter measurements were also done. The 
details are given below. 

2,3.1 Injector Flow .Analysis 

LABEOC VI operating manual provides the injector 
flow calibration curve for oxygen-methane propellant system. 

As shown in Appendix A.1? these can be deduced for oxygen - 
hydrogen propellant system. 
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FIG.9_ THRUST MEASURING UNIT 
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2.3.2 Flame-in limits Establishment 

Here for various fixed settings of either fuel or 
oxidiser, the setting of the other one is so varied until the 
flame just sets in. Combustion chamber pressure is measured at 
this instant. From regulator settings and injector calibration 
curves for hydrogen oxygen combination, 0/P ratio is calculated. 
Fixed oxygen setting gives high 0/F ratio flame-in limits while 
fixed hydrogen setting gives low 0/F ratio flame-in limits. 

This is done so as to obtain readings with maximum control para- 
meter provided with LABROC VI and to avoid dangerous explosions 
due to hydrogen leakage. .Appendix A. 2 deals with data., sample 
calculations and calculated readings. 

2.3.3 Thrust Meter Analysis 

Here, after the thrust knob is made to touch strain 
gauge plate with initial load, null point is adjusted on strain 
gauge indicator. Strain readings are taken for each incremen- 
tal load. Appendix A.3 deals with data, ssaple calculations and 
calculated readings. 

2.3.4 Determination of Performance Parameters 

For a particular mass flow rate through a nozzle, 
various regular readings for hydrogen and oxygen are fixed up 
for different 0/F ratios. The thrust meter is then calibrated 
as indicated in 2.3.3. After the engine rurfe, firing is made 
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according to instruction manual <.f LABR.0C VI. Readings of 
thrust meter, combustion chamber pressure, coolant flow rate 
and inlet and outlet temperature, etc. are taken at fixed 0/P 
ratio, for a given mass flow rate. The similar readings are 
taken at different 0/F ratios at the seme mass flow rate and 
at different mass flow rates. The same is repeated with diffe- 
rent nozzles. Appendix A.4 deals with data, sample calcula- 
tions and calculated readings. Tables ( 3 ), (4), ( 5 ) end (6) 
are used to determine the necessary Panel settings. 



CHAPTER III 


RE SUIiES AND DISCUSSION 

3.1 Injector Plow Analysis 

3.2 Plane-in LimrtB Establishment 

3.3 Thrust Meter Analysis 

3.4 Performance Parameter Analysis 
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3.1 Injector Plow Analysis 


A plot of the data, given in Table (l), Appendix A. 1 
is shown in Pig. (ll). This figure is obtained from the opera- 
tional manual of LABROC VI. Also a plot of the data given in 
Table (2), Appendix A.1 is shown in the Pig. (l2). 


Prom the Pig. ( 1 1 ) and Pig. (l2) it can be seen 
that when Panel pressure in Psig vs. flow rate in lb/ sec. is 
plotted on LOG x LOG scale, their variations are lineal'. This 
is in good agreement with the equation (A. 1 . 1 ) considered on 
LOG x LOG scale for incompressible flow 


m = C , A. / 2A/\p' 
d m /v.' ' t-** 


t \ 

or m = Const A P 


Taking logerithm on both sides we have 

. . log (m) = log const + l/2 ! ' log A p 

for p^_ /p Panel this can be written as 

. . log (m) = log const + l/2 log ( p Panel) 

Also for methane, molecular weight is 16 gms. /moles 

while for hydrogen it is only 2 gms. /moles. Therefore for the 

same panel setting less anount of hydrogen gas is obtained due 

to difference in molecular weights and hence densities, for 

fi and A. remaining the sane, 
d in 



PANEL PRESSURE 



FLOW RATE LBS/SEC. 

FIG. 1 1 - INJECTOR FLOW CALIBRATION FOR GASEOUS 
PROPELLANTS METHANE AND OXYGEN 
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3.2 Flame-in Limits Establishment 

A plot of the data given in the Tables (9) and (lo) 
in Appendix A.2 for a particular rocket motor, nozzle combina- 
tion is given in Fig. (13) using Hg - Og fuels and oxidiser. 

Also a plot of the data given in the Tables (13) & (l4) in 
Appendix A.2 for other rocket motor, nozzle combination is gi- 
ven in Fig. (l4) 

From the Fig. ( 13 ) and Fig. ( 14 ) it can be seen 
that a plot of combustion chamber pressure vs. 0/F ratio esta- 
blishes a zone for given rocket motor and nozzle combination. 

That is to say that the rocket motor must be operated within 
that zone to obtain the flame. It is also clear that this 
operating zone is quite wider for rocket motor and nozzle £ 1 1 1 
than for rocket motor and nozzle € 2.09 J 1. 

3.3 Thrust Meter Analysis 

The plots of datas given in Tables (15), (16), (l7), 
(iS ) , (19) and (20) in Appendix A. 3 are given separately in 
Figs. (15), (16), (17), (18), (19) and (20) respectively. 

From the Figures (15), ( 16 ) , (l7), (l8), (l9) aid 
(20) it is clear that strain developed in the load cell designed 
is linearly proportional to the load qpplied. This also indi- 
cates that all the time, we are working within elastic limit of 
the gauge. For the same initial load of 2.11 lbs., for diffe- 
rent sets, the initial end therefore the rest of the readings 



\> c KGF/CM 


ROCKET MOTOR 202253 SN/14 
5 — ROCKET NOZZLE 20226! SN/2! 
E 2.09'.l , D t =0.360CMS. 



O/F RATIO 

FtG.I 4 _FLAME-IN LIMITS OF GASEOUS PROPELLANT ROCKET 
MOTOR USING HYDROGEN AND OXYGEN 




STRAIN IN JU. IN /INCH 


4800 


4700 


4600 


4500 


4400 


4300 j 





4 


5 


LOAD IN LBS. 

FIG. 15-THRUST METER CALIBRATION CURVE 



STRAIN IN A IN /INCH 


5000r 


4900 


ROCKET MOTOR 2022E3 SN/14 

ROCKET NOZZLE 202726 SN/4 

L!'. I D t * 0.381 CMS 

(USED FOR PERFORMANCE - 
ANALYSIS SFT TWO) 


4800 


% 4700' 


4600 


4500: 


O I 2 3 4 

LOAD IN LBS. 

FIG I6_ THRUST METER CALIBRATION CURVE 


STRAIN IN M IN/INCH. 


33 



FIG.I7_ THRUST METER CALIBRATION CURVE 
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FIG. 19 THRUST METER CALIBRATION CURVE 
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are slight different. This is probably due to the following 
reasons s 

(1) Unequal temperature conditions 

(2) unequal tensions in the strings 

( 3 ) Accidental blasts which happened on account 
of leakage of hydrogen gas from the solenoid 
operated valve and 

( 4 ) Velocity of air vhich may shift the readings 
due to oscillations of hanging weights. 

3.4 Performance Parameter Analysis 

The plots of ihe data given in Tables (22), (24), 
(26), (28), (30) and ( 32 ) in Appendix A.4 are given in ihe 
Pigs. (21), (22), (23), (24), (25) and (26) respectively. 
These figures include the variation of performance parameters 
of gaseous propellant rocket motor using hydrogen and oxygen 
for 0/P ratio varying from 4 to 8 and for mass flow rates 
varying from 0*009 lb/sec, to 0*013 lb/ sec* 

Prom the Pigs. (2l), (22), (23), (24), (25) and 
for a particular set the following can be concluded . 

Ihen O/F ratio is decreasing 

(1) Combustion chamber pressure increases 

(2) Thrust and hence specific impulse also increases 

(3) Characteristic velocity, which is the measure , of combus' 
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(4) Equivalent velocity also increases. 

( 5 ) Thrust coefficient which is the measure of nozzle working, 

reduces. The reason being the '•count by which p builds' up 

c 

being more than build up .-cnuat n f. P. 

Conclusions (l) to ( 4 ) can be explained as under. Lower 
0/P ratio implies the higher settings of panel pressures. This 
leads to higher mixing pressure of hydrogen oxygen jet before 
Ignition. On firing obviously this leads to higher pressure 
which indicates the availability of more amount of energy for sub- 
sequent conversion. 

How we will observe in detail Tables (22), (24) and (26). 
It is clear that for increase in mass flow rate thrust increases 
(exception set 3 for its first two readings). But the amount by 
which, thrust increases is less than that of mass flow rate in- 
crease, resulting in the reduction of the specific. impulse (ob- 
serve last two readings of set 2 and for whole of set 3). Thrust 
coefficient variation range is largest for set 2 then comes for 
set 1 and set 3. Thus a kind of saturation occurs with set 3. 
Characteristic velocity for set 2 is larger than for set 1. Sane 
is -the case with then with equivalent velocity except last two 
readings. Equivalent velocity for set 3 is considerably smaller. 
This indicates that increase in thrust is not proportionate with 
mass flow rate increase. Higher characteristic velocity for set 
2 and set 3 is indicative of larger increase of chamber pressure 
than increase of mass flow rate. 
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On the basis of the higher expansion ratio nozzle used 
we can draw the following conclusions. Observe Tables (28), (30) 
end (32). 

The thiust F is increasing with increased mass flow 
rate. Specific impulse is improved for set 5 than for set 4 (ex- 
cept for reading 2), reason being thiust increases by larger 
amount than mass flow rate increase. But for set 6 once again the 
increase of mass flow rate dominates over the increase of thrust, 
resulting in lesser specific impulse than set 5. Thiust coeffi- 
cient variation range is maximum for set 5 then comes set 6 and 
set 4. Characteristic velocity aid equivalent velocity are higher 
for set 5 than for set 4 except for second reading in equivalent 
velocity oomparision. Characteristic velocity aid equivalent 
velocity are considerably smaller for set 6 than for set 5 (ex- 
cept for first reading), reason being the larger increase of mass 
flow rate than the increase of both chamber pressure end thrust. 

For both the nozzles optimum performance (both specific 

impulse and characteristic velocity) is obtained by operating the 

rocket motor at lower O/F ratios, of course at reduced nozzle per*- 

formance i.e.- thrust coefficient. For nozzle with 6=1*1 
2 

A . = TC / 4 D, is greater then for nozzle with £ - 2.09 * 1 
u "G 

A, by some 12$. The sane caa be the cause of reduced performance 

V 

with 2nd nozzle. For this compare specific impulse and thiust 
coefficient for set 1 and 4 (exception first two readings in spe- 
cific impulse comparison) and set 2 and 5. 
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1.2 Work Related To Theoretical" Performance 
Analysis 


1.3 Present Work; 
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1.1 General 

Theoretical investigations of performance analysis of 
chemical rockets are of prime importance to guide propellant 
research by leading the lines of efforts towards high perfo ro- 
mance system, to aid engine design by providing information on 
temperature, pressure, composition of product gases and optimum 
mixture ratio and to provide data input for flight performance 
calculations. All these prehand data can be made avail shle and 
a systematic theoretical approach to the problem of performance 
analysis is made. To handle the problem easily, appropriate 
assumptions are made. 

1.2 Work Related to Theoretical Performance Analysis 

In their report supported by HACA, Huff et al ( 1 9) 

S P 

attempted to device a general method for comutation of equili- 
brium composition and adiabatic temperature of chemical reaction 
with 3ST number of components in product gas composition. Beighley 
(20) and Donegan (21 ) also tried to use computer systems. 

Donegan tried to determine thermochemical propellant calculations 
on high speed digital computers while Beighley concentrated on 
reduction of rocket motor performance data using IBM computing 
machines. All their individual efforts were towards more general 
system of propellants. Morgan (22) applied theoretical approach 
for thrust determination of rocket motor for C-H-H-O— F system. 
Gordon et al' (23) worked car problem of theoretical performance 
analysis of chemical rocket motor using liquid oxygen and hydrogen. 



49 


In this technical note supported by NASA, Gordon (24) et al 
generalised methods for both equilibrium composition computa- 
tion end for theoretical rocket performance. They developed 
computer programme in *SQAP II', for the same. Moffatt et al (25) 
dealt with in particular H 9 - 0 ^ combination and determined pro- 
duct of combustion . at elevated temperature. Mentz (26) developed 
the computer programme for the subject concerning equilibrium 
temperature and composition determination. He nicely discussed 
the convergence techniques for the problem. Stanely (27) deve- 
loped an integrated programme for propellant performance analysis. 
All necessary data with graph required were' arranged within his 
programme. Olson (28 ) studied both theoretically end experimen- 
tally the recombination and condensation process in high area 
ratio nozzles for various propellents. Many standard, text books 
like Penner (29), Wilkins (30) have also dealt with this problem 
as a subject. Throughout Scarbrough's (31) book has been a use- 
ful reference for mathematics involved in the problem. 

1.3 Present Work 

The theoretical approach to the problem of performance 
analysis of a chemical rocket is initiated to have a method of 
computing from a given propellant combination of Hydrogen-Oxygen- 
Nitrogen system and chamber pressure; the adiabatic combustion 

temperature, equilibrium gas composition and other performance 

l 

parameter from both frozen flow and equilibrium flow techniques. 
Ultimate aim is to know the optimum working conditions. 



CHAPTER II 

THEORETICAL CAiOJLAIIOlTS 


2.1 Determination of Combustion Temperature and 
Composition of Product G-as 

2.2 Assumptions for Theoretical Performance 
Analysis of A Chemical Rocket and Introduction 
to Performance Parameters 

2.3 Frozen Flow Approach 

2.4 Equilibrium Flow Approach 


* • 
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2#1 Determination of Combustion Temperature and 
Composition of Product G-as 


Consider the following basic equation for chemical 
reaction of - 0 ^ - Ng system where is introduced as an 
impurity. 


A (m IT + o„ + K N„) n , H 0 + n 0 0 H + n, IT 0 

+ X4 H 2 + n 5 ° 2 + n 6 K 2 + n 7 H + n 8 0 + n 9 1T (2.1.l) 

where 

m = 16/OPE and 
K = (( 16/OPE) +l) x 0.001 

OPE = 0/P ratio. 

Aim is to determine A, (i = 1, 2, 9) and 

combustion temperature T. To simplify the calculations A and 
n^ are determined such that 

p.=n. (2.1.2) 

*x x 

i.e. species volume is specified to be ft T. Applying Dalton 1 s 
law of partial pressure we have for all gaseous phases 


P — + Pg — — + p^ — n^ t ng — — + n^ (2.1.3) 

law of conservation of mass is applied to the moles of species H 

a = -i- (2 n 1 + n 2 + 2 n 4 + /n ? ) = 2 m (2,1.4) 

for moles of species 0 

b = (n^ + ng + rij + 2 + n^) = 2 (2.1.5) 
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for modes of species N 


r ( v 2 vV = 2 K 


(2.1.6) 


Also consider the following dissociation equations 


h 2+ ^0 2 -h 2 0 


h„ 0 H 2 + 0 H 


^1 ( \l/2 

n 4 K 5 

n 2 (n 4 ) 


(2.1.7) 


(2.1.8) 




2 H 2?' H 


^ = (n 6 ) 1/2 (n 5 ) V ^ 


(2.1.9) 


(2.1.10) 


1 ° 2 f-° 


Rc = : 

5 (n 5 ) 


( 2 . 1 . 11 ) 


i b 2 E 


(2.1.12) 


Also .for adiabatic combustion 


(h!v 1 - te;) 0 ♦ K = I jSA 2 o n i + (H o >« a 2 + 


Thus we have a 


+ ^ E o\ n 9 j (2.1.13) 

svstem of eleven simultaneous nonlinear 


algebric equations from (2.1.3) to (2.1.13) for eleven untaomns 

to be determined, vis A, T end n ± (i = 1, 2, 9) for a given 

value of combustion pressure Old 0/F ratio. 
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-Analytical solution of eleven nonlinear simultaneous 
algebric equations being extremely tedeous, suitable technique 
may be used to solve them numerically. Newton - Rapson's method 
of iteration is used to make the equations linear which are then 
solved using matrix inversion method by Gauss - Jordan and next 
iteration is started if any. To speed up the iteration process 
natural logaxithm system is used. 


Consider equation (2.1.4). Take the natural logarithm 
of both the sides 


. . In a = In (2 n^ + n 2 + 2 n^ + n^) - In A 

Considering first order variations we have 

. 

h -‘ « 

/v win a a _ , - In a a - „ , 

3^ a - ^ La In n + - ln n A 1x1 n 2 + 

U 1 w 2 


"A 

Pin a 
p v ln n 

d 


Ain n. 


O' In a 
TMn n_ 


also 


i,J In a 


In n 


tv — t~ A In A 
\la. A 

J 

l wln a 


0 


1 


a n i 


d 


C.J n i 

p,ln n 1 


A In 


n„ + 


In a , 


d n i C 


• '“''•In n . 


cl n i 

2 n. 


= ( 


2 n 1 + n 2 + 2 n 4 + 


> n i = TTXJ 


(2.1.14) 


(2.1.15) 


Calc 


Proceeding similarly we have 


Oln a 

n 2 

O' In 

_a 

0in 

d Calc 

pita 

n 4 

a 

w In a 

n 7 . . 

01n 

_a 

d 111 d 

(a A) Calc 

0 ta 

A 


2 n 4 
(a a) 


Calc 


1 etc 
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also Ain. a = In a_ - In a_ _ 

IP rue Calc 

Thus Eg. (2.1,15) becomes 


/ ^xue \ 

= In ( — ) 

c 'Calc 


• • (a A >Calc 111 ( 


"Calc 


) = 2 A In n 1 + 


n 2 A In n 2 + 2 n 4 A In A In n^ 


- ^ A Wc A 111 A 


( 2 . 1 . 16 ) 


Similarly rest of the equations can be linearised. 
They become 


^ b A ^Calc 111 ^ 


Ture 

°Calc 


) = n 1 A In n 1 + n g A In n 2 + 


n^ Ain n^ + 2 n 5 A In n^ + n g A In n Q 

“ A ^Calc A l 11 A 


(2.1.17) 


A A ^Calc ln ^ c 


True 


Calc 


) = n_ A^n n + 2m A In n g + 


n Q Aln n 9 " ( ° A) Calc A 


(2. 1.18) 


P n l* 1 ( 

Calc 


'True 
' Calc 


i = 9 


) = X. n i 


A In n. 


(2.1.19) 


is 1 


On = Ain n - ilnn --f Aln« s (2.1.20) 

- “I — * 1 


In ( 


1 Calc 
®2 True 


E, 


2 Calc 


) = Ain n 2 + ~ Ain n - 1 . Ain n 1 


( 2 . 1 . 21 ) 


< ftx ) = Alnn S (2.1.22) 
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In ( 
In ( 

In ( 


B, m 

4 True 

R 4 Gale 


E. 


5 True 


E 


5 Calc 


E 


'6 True 


E, 


) = -- 


.Ain n^ + ,A In n^ 


) = - 4 A in n r + A In 


) = - i A in n, + /^ln n c 


6 Calc 
( h - 111 ( 


h, 


True 


Calc ~ ' h _ 

Calc 


-) = r, 

1=1 

i= 9 


where 


(h t ). 

o'x 


I 


i=1 p i 1 


C° dT + H° 
P O 


(2.1.23) 

(2.1.24) 

(2.1.25) 


(h A) 0alo A in A + S < I Ami (2.1.26) 


Note : P s= p . and T = T here* 
c c 

Thus modified equations from (2.1.16) to (2.1*26) are to 
be solved for Ain n 's, A In A and Ain T. Matrix form of 
the above equations is represented as on the next page. 
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*3 
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*8 
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1 

0 

0 

-1 
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0 

O O 
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1 

0 

*h 

0 

0 0 

O 

a 0 
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Ain *H 

— 


0 

o 

1 

0 

-Va 

-Va 0 

O 

0 0 

0 





0 

0 

o 

-Vi 

0 

0 1 

O 

0 0 

0 


A t>»hg 


^itrSto 

0 

0 

o 

0 

-V* 

0 0 

1 

a 0 

c 


AUt 9 


h. (s75 Tt.if \ 

x ^WcSi } 

o 

o 

0 

0 

o 

-V2. 0 

O 

1 0 

0 


AU a 




. »&a 


u \\ 

h 1 

>55 5 


ijyt 

H r -» C-^a) 

n ^9 Calf P w 


aUt ! 














_ — 


— — 


( 2 . 1 . 27 ) 


For numerical solution it is essential that diagonal elements 
should be non zero, necessary cure is taken in arranging the system 


as shown in Appendix B. 1 
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To start with the solution of the above system assume 
initial guess for n ±s A and T, Calculate all true thermo dynamic 
quantities at temperature T. Calculated values are found from 
the formulation. After finding the solutions, vhich are being 
error terms, correct the assumed values of n^, A and T, till 
the maximum error is some minimum quantity. Note that n.. = n^/A 
since value of A was so calculated to simplify the solution. 
Data, sample calculations, computer programme with necessary mo- 
difications for numerical solution and calculated data are re- 
presented in Appneix B. 1. 

2.2 Assumptions for Theoretical Performance Analysis 
of A Rocket Motor and Introduction to Performance 
Paranenters C See Pen-net- co.'p) 

Analysis is based on the following assumptions.. 

(1) Combustion process is adiabatic aid after it, thermodyna- 
mic equilibrium is assumed to be reached. 

(2) Products of combustion behave as ideal gases. 

(3) Expansion of products of combustion in a nozzle can be con- 
sidered as one dimensional flow of non viscous, ideal gases. 

(4) Expansion again is assumed to be adiabatic. 

(5) Velocity at the entrance of expansion nozzle is negligibly 
small in comparision to exit velocity from the nozzle. 

(6) Particular average value of y and W is assumed. 

With the above assumptions, following parameters for 

a chemical rocket motor can be defined. We will restrict to a 
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case of atmospheric expansion in a small divergence angle 
nozzle. 


1 v. 

Thrust P = ^ (l + cos a ) w — + (p, - p ) A 


a' e 


P = w u /g 
eq' ° 

vhere u is Equivalent or Effective exhaust velocity. 


u = ~ (l + cos CL ) v + 
eq 2 e 


K ~ Pa 1 J a 

w/g 


for p = p and Ot very small we have 

6 £L 


U = V 

eq e 


P = w v g /g 


Specific Impulse I = P/w = u eQ / g 


for p = p and 0< very small we have 
e a 


n -• A 

I = v /g v . 

sp e & / — 1 

— 

Also T /T = (p /pj ^ 


e c 


e' A c' 


„ g J RT 

Characteristic velocity C* = \ — * — — * ) 

_ r •* y 


where 


v + 


F* = ( irl — ) 

& y + 1 


2 (y*- i) 


Thrust Coefficient C^ = u g ^ / C* 


( 2 . 2 . 1 ) 

( 2 . 2 . 2 ) 

(2.2.3) 

(2.2.4) 

(2.2.5) 

( 2 . 2 . 6 ) 

(2.2.7) 

( 2 . 2 . 8 ) 

(2.2.9) 

( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 


P 


°F P = \ 


( 2 . 2 . 12 ) 
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Particulat average value of 


y 




y sad W is taken as 

— * 


7 


end 


W = W = W = W * 
c e 


(2.2.13) 


(2.2.14) 


2.3 Frozen Flow Approach 

Here recombination process occurs in a time i/hich is 
relatively larger in comparision to the expansion time. That 
is to say that reaction rates are very slow and the propellant 
products are said to be frozen at combustion chamber composition. 

/■pp lying law of conservation of energy at entry aid 
exit of rocket nozzle 


4 W v 2 = Ah® 
2 e c 


where 


j=9 

W = I, ¥, end X, 

5=1 3 3 


= n 


3 ^^Total 


also 


Ah® = T' s i 

3=1 3 3jI c 


w • ) 

3»®. 


(2.3.1) 


(2.3.2) 


(2.3.3) 


Now to determine T , we will assume the flow through 

6 

nozzle 'as isentropic flow though adiabatic flow will give the 
same results as developed by Penner (29). Thus for isentropic 
flow through nozzle we have 

S T C , P c = S ^ e > P e 


(2.3.4) 
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for mixture of ideal gases Penner (29) gives 

d=9 0=9 


T,P 




0=1 3 

Prom (2»3.4) and (2.3*5) we have 


In p + 


Y 


0=1 


X. In X. 
0 0 


(2.3.5) 


2? X, ( s° „ - s° ) = H m P /p 


j=1 3 3,T c 3,I e 


c A e 


(2.3.6) 


Brief procedure for performance parameter dete ruination is out- 
lined here. 


(1) Bind combustion temperature T = T c and composition using 
the theory developed in 2.1. 

(2) For p e = p Q = 1 at a, find T e by trial end error method 
using Eq. (2.3.6). 

(3) Find W using Eq. (2.3.2) 

(4) Find enthalpy change during expansion using Eq. (2.3.3) 

(5) Find v g using Eq. (2.3. l) for small CL and P e = P a 

(6) Find I using Eq. (2.2.7) for small 0 and P a = P„ 

sp © cl 

(7) Find y* using Eq. (2.2.8) 

(8) Find 0 * using Eq. (2.2.9) end (2.2, 10) 

(9) Find Cp using Eq. (2.2.11) 

Jill these calculations are for a fixed combustion cham- 
ber pressure p and for fixed 0/P ratio. The same is repeated 
c 

for various p 's and O/F ratios. Data, sample calculations, 
c 

computer programme aid calculated readings are given in 
.Appendix B.2. 



2.4 Equilibrium Plow Approach 
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Here reaction rates are too fast and various recombi- 
nation processes occur in a time which is short in comp arision 
to expansion time. Thus corresponding to each temperature thermo- 
chemical equilibrium is achieved and new composition determination 
is a must (or reaction rates for both forward and backward reac- 
tions are same at each temperature encountered in expansion flow). 


Equation (2.3. l) is still valid for equilibrium flow 
analysis except that one must restrict the analysis to a given 
weight of gas mixture, rather thaa to the temperature dependent 
molecular weight of gases. Choosing W q as Initial molecular wei- 
ght and then applying law of conservation of energy at nozzle 
entry aid exit we have 


i_9 W„ 

1 “ v 2 = H - H ( -2- ) 


2 "c e 


W 


e 


where 


W 


A;3 ’ V p c 

w = jp z. m w. 

e T e* P e J 


& 


(2.4.1) 


(2.4.2) 


(2.4..3) 


For H and H , absolute molar enthalpies of the gas mixture 

C 6 

referred to 298.16° K as reference point. Thus 
3=9 


h = y; x. m 

0 3=1 3 > *o p c 


“j + h ;, * 0 - v*.*] 


(2.4.4) 


and. 
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He " dll +H i, T e " H 3, 298.16 


(2.4.5) 


The ratio W in (2.4. l) corrects the molar enthalpy 
at nozzle exit to the enthalpy for the seme -/.’eight of gas mix- 
ture for which H is calculated, 
c 

For equilibrium flow the only method utilising the con- 
dition that the expansion is isentropic, is used for determina- 
tion of T g . Isentropic flow occurs for total weight of the gas 
and not for temperature dependent molecular weight of the mixture. 


s m = (i Ai ) s_ 

T , p c T ,p^ 

c 7 r G e 7 G 


(2.4.6) 


using Fq. (2.4.6) end Eq. (2.3.5) we have 


S X i T B ( S i T - * m p c - ® 111 X 3,t ,p ) 

An P„ 3» i c 0**13 


O' ^ o 




3,T - Slnp a - RtoX 3,T e ,P e ) 


(2.4.7) 


#iich can be further simplified to 


£ Jv 0 ,p 0 w J.*.- HlnI 3.Vp. ) " <5 ' /, * )x 


V e .P e (S 5,I ‘ » 10 -t. 5 


j> T e» p e 


R In, p^ - ( W n / W fi ) In p 


c e 


(2.4.8) 
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Equation (2.4.8) can be used to find*T by trial and 

e 


error method. Eor every T , p set; X. _ is to be determined. 

e’ *e 3,T o) p o 

— t/ t? 

Also y has no meaning in shifting equilibrium flow. Though 

— ■¥: 

it can be approximated by Eq. (2.2.8)- also V.' is modified to 

V* = j (W c + W e ) (2.4.9) 

Brief procedure is outlined here. 


(1) Bind T and composition using the theory developed 

c 

in 2.1 

(2) Eor p = p , calculate T by trial and error method. 

6 St 6 

Eor each assumed value of T find composition i.e. 

X. m 
^ T e’ p e 

( 3 ) Bind enthalpy change using Eqs. (2.4. l), (2.4.4) and 
(2.4.5) by 


W _ d=9 

H c -tr-)H e = £ 

W e & 


I . _ ( A Hf ° + H 

D,T c ,p c 3 


3,T c " H 3,298.16- 


W 

- ( — ) X . _ ( A Hf ° + H . 

“ ^ T e» p e 3 


W 


_ TT V 

3,T e n 5,298.l6 ;i 


(2.4.10) 


( 4 ) Bind v using Eq. (2.4. l) and (2.4. 10) 

( 5 ) Bor p = p and small Ci find I using Eq. (2.2.7) 

7 e a sp 

(6) Bind y* using Eq. (2.2.8) 

( 7 ) Bind C* using Eqs. (2,2.9), (2.2. 10) and (2.4.9) 

(8) Bind using Eq. (2.2.11) 

All these calculations are for a fixed combustion chamber 
pressure and 0/B ratio. The sane is repeated for various p^s and 
0/E ratios. Data, sample calculations, computer programme and 
calculated data is given in Appendix B.3. 


CHAPTER III 


HESUHES AID DISCUSSION 


3.1 .Analysis of Combustion Temperature and 
Composition of Product Gas 

3.2 Frozen Flow Approach 

3.3 Equilibrium Flow Approach 


5» 1 Analysis of Combustion Temperature and 

Composition of Product &as 

The plots of the data given in the Tables (53), (34), 
(35) and (36); (37), (38), (39) and (40); (4 1 ) , (42), (43) and 
(44); (45), (46), (47) and (48) are given in Figures (27), (28), 
(29) and (30). These figures show the variations of combustion 
temperature and mole fractions of various species assumed to be 
present in the product gas, for combustion chamber pressures 
varying from 50 Psia to 200 Psia, in step of 50 Psia and O/P 
ratio 2 to 8, in step of 2 each. 

Prom -the figures (27), (28), (29) and (30), for a 
particular set following can be concluded when 0/P ratio is in- 
creasing. 

(1) Equilibrium combustion temperature increases. At hi- 
gher 0/F ratio this increase is flatten. 

(2) Mole fraction for water increases. At higher 0/P ratio 
this increment is reduced. 

(3) Mole fraction for OH increases. 

(4) Mole fraction for 10 increases. 

(5) Mole fraction for H 2 decreases. 

(6) Mole fraction for 0 2 increases. 

(7) Mole fraction for N 2 decreases. 

(e) Mole fraction for H first increases and then decreases. 

(9) Mole fraction for 0 increases. 

(lO) Mole fraction for H first increases and then decreases. 
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Conclusions which include the increase of quantity can 
be explained as under. Availability of more oxygen with increased 
0/3? ratio tends to complete combustion resulting in increased quan- 
tities of T q and mole fractions of HgO, OH, HO, 0 , 0 etc. Mole 
fraction of II ^ decreases with 0/3? ratio increase because lesser 
amount of H^ is available for taking part in reaction. Mole frac- 
tion of decreases as it is accounted for in factor K, in formu- 
lation which decreases with 0/3? ratio. Thus lesser available Ng, 
results in lesser N 2 tn product gas. Dissociated mole fractions 
of N and H increases with 0/3? ratio at lower values of 0/P ratios 
but at higher values of O/F ratio they decrease with increase in 
0/P ratios. Indicative of relative supression of dissociation of 
N^and H a at higher 0/P ratio. Higher value of combustion tempera- 
ture is due to the inclusion of all basic nine species, in parti- 
cular NO, 0 and N, in formulation of the problem. 

Prom the Table group (33), (37), (4l) end (45); (34), 
( 38 ), (42) end (46); (35), (39), (43), aid (47); ( 36 ), (40), (44) 
and (48) for increasing combustion chamber pressure, at particular 
0/P ratio, following cm be- concluded. 

(1) Combustion temperature .increases. 

(2) Mole fraction for H 2 0 decreases. 

( 3 ) Mole fraction for OH .increases. 

( 4 ) Mole fraction for NO increases. 

( 5 ) Mole fraction for B 2 increases. 

( 6 ) Mole fraction f or N 2 decreases. 
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( 7 ) Mole fraction for H decreases. 

(8) Mole fraction for N increases. 

(9) Variations in mole fractions of 0 2 end 0 are bit non 
predictable. 

Sole reason being the increase in combustion tempera- 
ture with 0/F ratio. 

3«2 Frozen Flow Approach 

The plots of data given in the Tables (49), (50), (’5l) 
and (52) are given in the figure (3l), (32), (33) aid (34). These 
figures include the variations of performance parameters of che- 
mical rocket with propellant system for combustion cham- 

ber pressures 50, 100, 150 and 200 Psia and 0/F ratios 2, 4, 6 and 
8 at each chamber pressure. 

From the figures (31 ), (32), (33) aid (34), for a par- 
ticular set, following can be concluded. TShen 0/F ratio is 
increasing 

(1) Exit temperature increases. 

( 2 ) Average molecular weight increases. 

( 3 ) Exit velocity decreases. 

( 4 ) Equivalent velocity decreases. 

(5) Specific impulse decreases. 

(6) Trend for y* aid thrust coefficient being little hit 
non predictable ♦ 

( 7 ) Characteristic velocity decreases. ; 



4000 6000 8000 









CHAMBER PRESSURE = 10.20 ATA 

= 150 PSIA . 
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The conclusions can be explained as under. With increa- 
sed 0/f ratio combustion temperature increases and nozzle being 
unchanged, exit temperature increases. Average molecular weight 
increases with increased and adjusted mole fractions of various 
species. Combined effect of above two reasons may be argued 
as the cause of conclusions (3), (4), (5) and (7). 

from the tables (49), (50), (51 ) and (52) following 
can also be concluded, for increasing combustion chamber pressure 
and at fixed 0/f ratio. 

(1) Exit temperature decreases. 

(2) Aversge molecular weight does not get affected a lot 
or it is approximately constant. 

(3) Both linear exit velocity and equivalent velocity 
increases. 

(4) Specific impulse increases. 

(5) ■y * variations are bit less predictable but it gene- 

rally decreases.. 

(6) Characteristic velocity increases. ' 

(7) Thrust coefficient increases. 

Thus it is obviously advantageous to run the rocket 
motor at fixed 0/f ratio with increased p c to obtain more effi- 
ciency of both combustion as well as nozzle performance. 

Also from Table (49) Sid figure <3l) it is clear that 
at lower chamber pressure, both comhrstion process as well as 
' nozzle performance are not efficient. This is obvious and is ' 



indicated by lower values of both thrust coefficient end diffe- 
rence between equivalent velocity and characteristic velocity. 
Higher values of exit temperature end lower values of specific 
impulse for this case are also indicative of the same. 

3*3 Equilibrium Plow Approach 

The plots of data given in the Tables (53), (54), 

(55) and (56) are given in the Figures (?5), (36), (37) and 
(38). These figures include the variations of performance 
parameters of a chemical rochet motor, theoretically calcula- 
ted for Hydrogen-Oxygen-Hit rogen propellant system, with equi- 
librium flow approach; for combustion chamber pressures of 
50, 100, 150 and 200 Psia and 0/P ratios 2, 4, 6 saad 8 at each 
chamber pressure. 

Conclusions drawn from these figures are more or less 
similar to those mentioned in Section 3*2, Frozen Plow Approach. 
IHhile comparing the results obtained from equilibrium flow 
approach with those obtained from the frozen flow approach, fo- 
llowing can be concluded. 

(t) At higher chamber pressures, performance parameters 
viz. specific impulse, characteristic velocity, equivalent ex- 
haust help city, linear exit velocity, average molecular weight, 
exit temperature etc* are having larger values than those -theo- 
retically calculated with the frozen flow approach. IhUe bolh 

■y * and thrust coefficient are bit lower. 



2000 4000 6000 0 M/ShC 
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(2) At lower chamber pressures, the theoretically calcula- 
ted performance parameters like specific impulse and both equi- 
valent exhaust velocity end linear exit velocity are lower than 
those by frozen flow approach, laile exit temperature, average 
molecular weight, characteristic velocity etc. are larger for 

. as. 

equilibrium flow approach case. Both. <y aid. thrust coefficient 
are lower for equilibrium flow approach case. For very lower 
chamber pressure of 50 Psia, thrust coefficient has reduced to 
even less than. 1. 

Reason for conclusion No. (l) is the additional energy- 
release, arises primarily from the recombination of dissociated 
and other species generated in combustion' chamber vhich are not 
stable at lower temperatures. This is prominent with higher 
chamber pressures. While at lower chamber pressures because of 
favoured dissociation of species like H^O, Hgj 0^, Hg P ar ” 
ticular temperature) this additional energy release may be redu- 
ced and hence reducing the total heat drop in nozzle. Moreover 
because the average molecular weight also increases the ratio 

Ah / W reduces, reducing both linear exit velocity and 
c 

equivalent exhaust velocity and hence specific impulse is also 
reduced. Increased characteristic velocities are due to lowered 
<y * 9 which overpowers molecular weight increase. 

ilso, for the case of lower chamber pressure, reduction 
of thrust coefficient is (tie to both reduction of equivalent 
exhaust velocity and increase of da arat eristic, velocity. This 



is prominent at chamber pressure of 50 Psia. Also note that 
— * 

V has actually no meaning with changing composition flow 
i*e« equilibrium flow because this contradicts, one of the 
basic assumptions. 

Prom Tables (53), (54), (55) sac! (56), following 
can be concluded as exceptional case, while making compari- 
sions with frozen flow conclusions. For increasing combus- 
tion chamber pressure and at fixed 0/F ratio, the average 
molecular weight here increases a bit. 



heoommehdadions 


Experimental set-up can be extended to study 

(1) Pressure distribution across the rocket nozzle. 

(2) Heat transfer characteristics of water cooled rocket 
motor. 

(3) Socket flame spectroscopy and spectrometry. 

Theoretical calculations can be extended with inclu- 
sion of condensed phases in product gas composition. Though 
theoretical propellant evaluation represents max imum per- 
formance and not necessarily that -which can be achieved in 
practice, matching can be tried with different compositions 
of product gas by suppressing some of the species for each 
case. 
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APPENDIX A 


DADA, SAMPLE CALCULATIONS AND PINAL 
CALCULATED HEADINGS POE EXPERIMENTAL 
APPROACH 


A. 1 Injector Plow Analysis 

A.2 Plane-in Limits Establishment 

A. 3 Thrust Meter Analysis 

A.4 Determination of Performance Parameters 
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A. 1 Injector Plow .Analysis 

Operational manual for LABEOC 71 provides injector 
flow calibration curve for SKT/14 unit. This includes plot of 
flow rate in lb/sec. Vs. panel pressure in Psig for oxygen 
and methane gas, on which it was originally designed. The 
following data is given by Pig, (ll). 

TABLE 1 

Injector plow Calibration Data for Gaseous 

Propellant Methane. 


GAS USED 

JPEMEL 

jPEESSUBE 
1UT Psig 

{PLOW RATE 

f nr 

Jibs/ S6 C « 


82 

0.00142 

Methane CH, 

4- 

170 

0.00280 

With Molecular 

Weight 16 gms/mole 

270 

0.00430 


The mass flow rate of the propellant sad the pressure 
drop across the injector are related by the equation 


» - °d P A in J ~ 2 A s/ p’ 


(A. 1.1) 


Por fixed C d , A^ and A P 

A = Const J]p = Const Ca.1.2) 

■ Por Hydrogen and Methane gas because we are selecting sane pre- 

, * • ,_1 . \ 

ssures and and hence \ / Pressure Head ' 




*1 / *2 



• * 


(A. 1.3) 
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Subscript 1 is for methane aid 2 is for hydrogen. Also W = 16 
gms/mole ahd W 2 = 2 gms/mole, 

’ * =a/~ 8 ' m 2 = 2.82843 m 2 

Y * 2 ‘ 8 3& 2 (A. 1.4) 

Thus at the same pressure settings we have 

(® 2 ) 1 = 0.00142/2.83 = 0.000502 

(m 2 ) 2 = 0.0028/2.83 = 0.000989 

(m 2 ) 3 = 0.0043/2.83 = 0.001519 

Thus we have following data for injector flow calibration for 
Hydrogen gas. 

TABU 2 

Injector Plow Calibration Data? Calculated 
For Gaseous Propellant Hydrogen. 



I PALED 

| PLOW RATE 

GAS USED 

? PRESSURE' 

l IF 


i IF Psig l Lb/Sec. 

Hydrogen H 2 

82 

0.000502 

With Molecular. 

170 

0.000989 

Weight 2 Gms/Mole 

270 

0.001519 


Thus similar graph for H 2 calibration caa be plotted. 
Same graph for 0 2 is used see Pig. (12). 

How problem is to find out panel pressure settings for 
both H g and 0 2 for various 0/P ratio aid at various mass flow 
rates through nozzle. 
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Take m = 0.009 lb/sec, through nozzle and 0/F ratio 

of 3 is selected. Row 

Rate of 0 2 Plow + Hate of H 2 Plow = m (A.1.5) 

» 

• • 4 (Rate of H 2 Plow) = 0.009 lb/sec. 

. . Rate of H 2 Plow = 0.009/4 = 0.00225 aid 
Rate of 0 2 Plow . . 3 x 0.00225 i.e. 

= 0.00675 both in lb/sec. 

Prom Pig. (12) panel setting forHg gas is 205 Psig 
and for 0 2 gas is 420 Psig. Same c cm bo repeated for various 
0/p ratio and at various m. Thus we have the following tables 
for calculated readings. 

TAKES 3 

Calculated Panel Pressure Settings for Caseous 

Propellants Hydrogen and Oxygen at Different 

0/P ratios form = 0.009 Ib/sec, 


0/P RATIO 

| PANEL j 

| PRESSURE j 

| IN Psig 

f I0B °2 j 

[ PANE! 

( PRESSURE 
( IN Psig 
| FOR H 2 

3 

205 

420 

4 

220 

350 

5 

230 

270 

6 

237 

225 

7 

242 

195 


8 


246 


170 



TABLE 6 
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Calculated Panel Pressure Settings For Gaseous 
Propellants Hydrogen and Oxygen at Different 
0/P Ratios form = 0.013 Ib/sec. 


0/P RATIO 

"liVl.n'L 
^PRESSURE 
IIN Psig 
|I0R 0 2 

| PANEL 
$ PRESSCIEE 
| IN Psig 
| FOR H 2 

3 

315 

560 

4 

328 

495 

5 

342 

405 

6 

350 

345 

7 

357 

295 

8 

395 

283 


Note s Actual sets of readings are restricted die to 

maximum range on regular being 600 Psig and due 
to leakage failure at high pressure regulator sett- 
ing. Thus 0/P ratio range is restricted. 

A.2 Flame-in Limits Establishment 

For rocket motor 202253 SN/14; Rocket Nozzle 
202726, SN/4, € 1 : 1, 3> t = 


Ctas 



TABLE 7 


Observed Readings Bor Elame-in Limits Establishment 
Experiment Eor Eixed Oxygen Setting. 


EDIT DO 




4 


EIXED SETT I MJ BOR 0 
IN' Psig 

0 2 SETTING AT EIRE, 
IN Psig 

H 2 SETTING- AT EIRE, 
IN Psig 

COMBUSTION CHAMBER 
PRESSURE AT EIRE, 


too 

200 

300 

400 

72 

163 

255 

340 

30 

43 

50 

60 

00 

05 

10 

16 


IN Psig 


TABLE 8 

Observed Readings Bor Blame— in Limits Establishment 
Experiment Bor Eixed Hydrogen Setting. 


RUN NO 


FT 

JLl 




EIXED SETTING BOR H 

IN Psig d 

100 

200 

300 

400 

H 2 SETTING AT EIRE, 

90 

185 

280 

375 

IN Psig 





0 2 SETTING AT EIRE, 

IN Psig 

10 

20 

30 

45 

COMBUSTION CHAMBER 
PRESSURE AT EIRE, 

IN Psig. 

00 

00 

00 

00 



Sample calculations are carried out for readings in 
Table 7. Prom Pig. (12) at pressure's noted fbrHg and 
their flow rate is calculated. 

At 72 Psig for 0^, flow rate is 0.00255 lb/sec. 

At 30 Psig for H , flow rate is 0.00019 lb/sec 

. . 0/P ratio = 0.00255/0.00019 

= 13.40 also 

2 

Combustion chanber pressure is 14.7 Psia i.e. 1.033 Kgf/cm . 
Calculated readings for data in Table (7) and Table (8) are as 
under. 


TABLE 9 

Calculated Data Por Plane-in Limits Establishment 
Experiment Por Pixed Oxygen Setting 


RUN NO 

1 : i : 1 

i 1 i ;i 1 

3 

T 

JL 

4 

HIGH 0/P RADIO 

13.40 20.00 

25.20 


28.80 

COMBUSTION CHAMBER 
PRESSURE AD EIRE, 

IN Fgf/Cm 2 

1.033 1.34 

1.68 


2.09 


TABLE 10 




Calculated Lata Por Plame-in Limits Establishment 
Experiment Por Pixed Hydrogen Setting 

HJN NO 

f 1 « 2 1 

3 

T 

JL 

4 

LOW 0/P RADIO 

0.788 0.755 

0.740 


0.730 

COMBUSTION CHAMBER 
PRESSURE AD EIRE, 

IN Kgf/Cm 2 

1.033 1.033 

1.033 


1.033 
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For rocket motor 202253 SN/lA, Socket nozzle 202261, 
SH/21, 6 2.09 s 1, D t = 0.360 Oms. 

TABLE 1 1 


Observed Headings for Flame-in limits Establishment Experiment 
For Fixed Oxygen Setting 


HUH HO 


I 

I 




4 


FIXED SETTING FOE 0 o 

IN Psig 

100 

200 

300 

400 

0 2 SETTING AT FIEE, 

IN Psig 

80 

175 

265 

360 

H 2 SETTING AT FIEE, 

IN Psig 

30 

100 

160 

225 

COMBUSTION CHAMBER 
PEESSUHE AT FIEE 

IN Psig 

00 

00 

AO 

62 


TABIE 12 

Observed Headings for Flame-in limits Establishment Experiment 
For Fixed Hydrogen Setting. 


1 ' 1 1 ,, J " 1 - — ■ 

iui: :to j 

r~r^ 

i 1 i 

m 

( i 

3 | 

4 

FIXED SETTING FOE H 

IN Psig 

100 

200 

300 

400 

H 2 SETTING AT FIEE, 

IN Psig 

85 

175 

26 5 

350 

0 2 SETTING AT FIEE, 

10 

25 

50 

70 


IN Psig 

COMBUSTION CHAMBEE 00 00 00 

PEESSUHE iff FIEE, 

IN Psig 


00 











Calculated readings for data in Table (ll) and 
Table (12) are as under. 


O' ABIE 13 

Calculated Data For Plane -in limits Establishment 
Experiment For Fixed Cxygen Setting, 


EDIT DO 

l , 1 
r 1 x 

2 -T 
2 11 

, nr 
3 1 

A 

*+ 

HIGH 0/F EAT 10 

14.75 

9.675 

8.825 

8.59 

COMBUSTION GE AMBES 
PEESSDEE ID Kgf/Cm 2 

1.033 

1.033 

3.72 

5.225 


TABLE 14 




Calculated Data For Plane-in limited Establishment 
Experiment For Pieced Hydrogen Setting. 

EUD NO 

1 1 f 

2 l 

3 ! 

4 

LOW 0/F RATIO 

0.83 

0.956 

1.30 

1.317 

COMBUSTION CH /IEEE 
PEESSDEE AT FIEE, 

ID Kgf/Cm 2 

1.033 

1.033 

1.033 

1.033 


A. 3 Thrust hoter Analysis 

From the set-up as shovoi in Fig, 10 it is "clear 
that how one should calibrate the thrust meter before taking 
set of performance parameters measurement. 
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For rocket motor 202253 } Si!/l4; Locket nozzle 202726 
SIT/4 , f 1 ! 1, D t = 0.381 Oms. 

TABLE 15 

Observed leadings For Thrust litter Calibration Experiment 
For I erformaice imalysis Set 1. 


LOADS IH 
Lbs 

l T , rT 

y k'JL . -j -Lj.< : 

Min/ In + ve 

2.11 

4400 

3.11 

4526 

4.11 

4676 

5.22 

4830 


TABLE 16 

Observed Headings For Thrust Meter Calibration Experiment 
For Ferformance .Analysis Set 2. 


LOADS IN I STEAD! Ill 

Lbs j U In/In + ve 

2.11 4554 

3.11 4696 

4.11 4832 


5.22 


4970 
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TABLE 17 


Observed Headings For Thrust Meter Calibration Experiment 
For Performance Analysis Set 3. 


LGAJD 

Lbs 

IN 

5 STRAIN IN 
| U- In/In + ve 


2.11 

4300 


3.11 

4628 


4.11 

4766 


5.22 

4908 


For rocket motor 2022~3j SN/14; Pocket nozzle 202261, 
SN/21, 6 2.09 : 1, D+ 0.360 Qns. 

TABLE 18 

Observed Readings For Thrust Meter Calibration Experiment 
For Performance Analysis Set 4. 


LC/D III 5 STRAIN IN 

Lbs | Mln/ln + ve 


2.11 

4740 

3.11 

4870 

4.11 

4992 

5.22 

5114 



TiiBEB 19 


Observed Headings For Thrust Meter Calibration Experiment 
For Performance .Analysis Set 5. 


1C id) III l STlilN II 
Lbs, \ U- In/ln -f ye 

2.11 4-760 


3.11 4877 

4.11 4975 

5.22 5094 


TABLE 20 

Observed Headings For Thrust Meter Calibration Experiment 
For Performance .Analysis Set 5. 


LOAD IN | ST BAIN IN 
Lbs. f Min/ In + ve 


2.11 

4770 

3.11 

4864 

4.11 

4982 

5.22 

5088 


ill thrust meter calibration readings are taken with 


cooling water circulation on. 



Am^'r De te min ation of F e rf ormance Parameters. 


For rocket motor 202253, SN/14; Socket nozzle 202726 
SU/4, 6 1 : 1, D t - 0.381 Oms. 

TABLE 21 

Co sewed Readings Ror Experiment of Determination of 
Performance Parameters of A Gaseous Socket Motor Using 
Hydrogen and Oxygen At Different 0/P Ratios For 
m = 0.009 Ib/sec. 

Performance .Analysis Set 1 


0/P RATIO 

}.P c Psig 

~~[lMR* 
ijU In/In 

H + ve 

(COOLANT 
(FLOW RATE 
Ilbs/Min, 

(COOLANT | COOLANT 

(inlet ( outlet 

(TEMP. °F I TEMP. °F 

* 8 

32 

4560 

15.20 

68 

84 

7 

35 

4570 

15.00 

68 

84 

6 

40 

4588 

15.40 

68 

88 

5 

46 

4600 

15.20 

68 

96 

4 

58 

4620 

15.40 

68 

108 


* Thrust meter readings with initial load = 2.11 lbs. 


For tliis case sample calculations are as under. 

O/P ratio is selected as 8 for sample calculation and 
all other data corresponding to it is taken. 

p 

Combustion chamber pressure in Kgf/CcT 
= (Observed reading + 14.7) / 14.7 
= (32 + 14.7) / 14.7 = 3.18 Kgf/Cm^ 
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Kow thiust calculations are to be done. 

Thiust meter reading is 4560 fJL In/In + ve. See Pig. (15) and 
read correspondingly the load of 3.32 lbs. Actual additional 
force = Thiust is 3.32 - 2.11 = 1.21 lbs. 

Thiust P = 1.21/2.2 = 0.551 Igf. 

Specific Impulse = P/m = = 134.5 Kgf Seo./Kgm 

Thiust Coefficient j/p 0 A t - vJvH ' olltt 0E 

_ --C .5 51 

0.114 x 3.18 “ u 

T^i e re A^ = 7T/4 ]y_ vhere 1^ = 0.381 Gms. 


Characteristic Velocity = p A + / m 

C *0 

32.2 x 46,7 x 0.114 
0.009 x 6.42 

= 3000 Pt/Sec. OR 

3.18 x 0.114 x 9.81 
0.009 / 2.2 

= 915 M/Sec. 


Equivalent Exhaust Velocity = 


Specific Impulse x g 
134.5 x 9.81 
1322 M/Sec. 


To get more correct results, carry out all calculations 
in EPS System and finally convert them in MKS System. , Though 
both systems are used here. Por data in Table (2l) calculated 
readings are tabulated as under. 



TABLE 22 


Calculated Lata Lor Experiment of Determination 
Cf Performance Parameters Of A Gaseous Pocket Motor 
Using Hydrogen and Oxygen s At Different 0/P Ratios 
Por r* = 0.009 Lt/Sec. 


o/p 

RATIO 

t: 1 

|Kgf/Cm ] 

FT" 

| Kgf. 

Ti ] 

f 

ftKgf Sec/ 
/Cgm 1 

! 

f 

T“ 

| M/Sec. 

T ~ — 

Q eq 

I M/Sec, 

8 

3.18 

0.551 

134.50 

1.48 

915 

1322 

7 

• 

v >3 

00 

0.586 

143.50 

1.465 

970 

1415 

6 

3.72 

0.636 

155.70 

1.440 

1060 

1540 

5 

4.13 

0.678 

165.50 

1.380 

1182 

1632 

4 

4.94 

0.745 

182.10 

1.270 

1420 

1800 


TABLE 23 

Observed Readings Por Experiment of Determination of 
Performance Parameters of A Gaseous Rocket Motor Using 
Hydrogen and Oxygen At Different Q/P Ratios Por 
m = 0.011 Lb/Sec. 

Performance Analysis Set 2 


0/P RATIO 

h 

1 c 
1 

Psig 

I TMR* 

| fJLln/ln 
| + ve 

.(COOLANT l OOOLJfT 
{PLOW RATI* INLET 
iLbs/Min. f TEMP °P 

( OOOLAN^ 
f OUTLET 
f TEMP. °P 

8 


50 

4796 

15.2 

64 

82 

7 


56 

4798 

15.2 

64 

84 

6 


61 

4802 

15.2 

64 

86 

5 


69 

4805 

15.3 

64 

90 

4 


79 

4806 

15.2 

64 

98 


*Thrust meter readings with, initial load = 2.11 Lbs. 

Use Pig. (l6) for calculations. Calculated readings 
for the ' data in Table (23) are as under. 



TABIE 24 


Calculated Data Dor Experiment of Determination And 
Of Performance Parameters Of A Gaseous Socket Motor 
Using Hydrogen and Oxygen, At Different 0/P Ratios 
Pbr m = 0.C11 Db/Sec. 


0/P 

RATIO 

i P 

* c . 2 
pgt/Cif 

T 

f 

if 

A 

P 

Kgf. 


JTi 

U P 

1 c*~ 

fM/Sec. 

I __ 

I u 

| M/Sec. 

8 

4.49 


0.792 

158.1 

1.51 


1035 

1551 

7 

4.92 


0.800 

160.0 

1.40 


1102 

1570 

6 

5.27 


0.809 

161.75 

1.33 


1212 

1585 

5 

5.82 


0.816 

163.00 

1.21 


1338 

1594 

4 

6.52 


0.821 

164.00 

1.05 


1500 

1609 


TABLE 25 

Observed Readings Por Experiment of Determination of 
Performance Parameters of A Gaseous Rocket Motor Using 
Hydrogen aid Oxygen At Different 0/P Ratios Por 
m = 0.013 lb/Sec» 

Performance Analysis Set 3 


0/P RATIO 

|p c Paig 

Jtmr* 

Xlllxi/ln 
$ + ve 

■ |COOLANT 
|pi£)W RATE 
libs/Min. 

JOOOlijffT 

lliCTT 

$TEMP °P 

jtCOOLAMP 

foUHET 

PIMP. °P 

8 

72 

4727 

15.2 

64 

80 

7 

73 

4730 

15.2 

64 

84 

6 

78 

4740 

15.2 

64 

68 

5 

87 

4744 

15.2 

64 

'■ 9& 


* Thrust meter readings with initial load = 2*11 Iifcte* 

Use Pig. (l7) for calculations; Calculated readings 
for data in Table ( 13 ) are as under. 
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Calculated Data Dor Experiment cf Determination .tod 
Of Performance Parameters Of A Gaseous Pocket Motor 
Using Hydrogen end Oxygen, At Different O/P Ratios 
Porm = 0.013/Lb/Sec. 


o/p 

RATIO 

} P c * 

{Kgf/Cm 2 I 

1 \ 

P ) 

Kgf. j 

pI7T 

Kgf.Sec/J 
f K£m 1 

C ' 

P < 

1 

< 

1 

[ c*~“ 

( M/Sec. 

( 

— " 

v eq 

0 M/Sec. 

8 

5.89 

0.782 

132.10 

1.120 

1172 

1299 

7 

5.96 

0.795 

134.50 

1.117 

1185 

1320 

6 

6.31 

0.836 

141,75 

1.115 

1253 

1390 

5 

6.94 

0.858 

145.50 

1.048 

1372 

1430 


Por rocket motor 202253, SN/14; Rocket nozzle 202261, 
SK/21, € 2.09 * 1, D. = 0.360 Cms 

TABLE 27 

Observed Readings Por Experiment of Determination of 
Performance Parameters of A Gaseous Rocket Motor Using 
Hydrogen aid Oxygen At Different C/P Ratios Por 
m = 0.009 Lb/Sec. 

Performance .Analysis Set 4 


O/P RATIO Jp 

Psig 

ItMR* fcOOLMT jfcOOLAUT 

I j U In/In fpLOW RATE^HUBT 

I + ve flbs/Min. $TEMP ®P 

I 000LAHT 

I OUTLET 

I TEMP. °P 

8 

40 

4894 

14.6 

66 

82 

7 

43 

4906 

14.5 

66 

86 

6 

48 

4910 

14.6 

66 

88 

5 

56 

4916 

14.8 

66 

96 

4 

66 

4930 

14.8 

66 

104 

* Thrust meter readings with 

initial lo ad = 2.11 

Lbs 


Use Pig, ( 18) for calculations. Calculated readings for data 
in Table (27) are as under. 



TABLE 28 


Calculated Data For Experiment of Determination And 
Of Performance Parameters Of A Gaseous Pocket Motor 
Using Hydrogen and Oxygen, At Different 0/E Ratios 
For el = 0.009 Lb/Sec. 


0/F 

RATIO 

T p c | 

kgf/Gm 2 | 

r» ] 
| K ^ f * | 

fl 
r ap 

Egf.Sec/ 

1 Egm j 

1 °* ■ 
i : 

r~c* ; 

| IS/ Se c . 

f e 9 

j K/Sec. 

8 

3.72 

0.5515 

136 

1.40 

947 

1335 

7 

3.92 

0.605 

148 

1.45 

1000 

1452 

6 

4.27 

0.623 

152.5 

1.37 

1086 

1495 

5 

4.82 

0.636 

156 

1.245 

1225 

1530 

4 

5.50 

0.723 

168 

1.176 

1400 

1640 


TABLE 29 

Observed Readings For Experiment Of Determination of 
Perfommoe Parameters of A Gaseous Rocket Motor Using 
Hydrogen and Oxygen At Different 0/F Ratios For 
m = 0.011 Lb/Sec. 


Performance Analysis Set 5 


0/F RATIO Jp Psig 

:__r 

l tmr* 

| jXln/ln 
i + ve 

icOOLAST J COOLANT 
iFLOW RATE* INLET 
Uibs/lin . t TEMP °F 

l COOLANT 
| OUTLET 
f TEHP. °F 

8 

54 

4928 

15 

67 

84 

7 

65 

4934 

15 

67 

90 

6 

74 

4945 

15.2 

67 

92 

5 

80 

4952 

15.2 

67 

94 


*Th.rust meter readings with, initial load = 2.11 Lbs. 

Use Figure (l9) for calculations. Calculated readings 
for data in Table (29) are as under. 
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I ABIE 30 

Calculated Data Dor Experiment of Determination And 
Of Performance Parameters Of A Gaseous Rocket Motor 
Using Hydrogen aid Oxygen, At Different 0/P Ratios 
Form = 0.011 Lb/ Sec. 


0/F 

RATIO 

— f 

\ Kgf/Cffi 2 f 

_L_ I 

F 

Kgf. 

n — r 

Y SP | 

| Kgf .Sec J 
1 Kgm i 

°F 

1 c*~ 

| M/Sec 

l u 
jt eq 

* l M/Sec . 

I 

8 

4.68 

0.708 

141.9 

1.425 

976 

1390 

7 

5.43 

0.731 

146.2 

1.250 

1132 

1432 

6 

5.96 

0.776 

156.3 

1.215 

1262 

1532 

5 

6.45 

0.804 

161 

1.175 

1349 

1585 


TABLE 31 

Observed Readings For Experiment of Determination of 
Performance Parameters of A Gaseous Rocket Motor Using 
Hydrogen and Oxygen At Different O/F Ratios For 
m = 0.012 Lb/Sec. 

Performance Analysis Set 6 


O/F RATIO I 

p Psig 
c 

I TIE* • 

1 /Xln/In 
i -fve 

I COOLAHT 
| FLOW RATE 
*f Lbs/Kin . 

l COOLANT 
| INLET 
| T1MP °3? 

J 000LJTT 
| CUTLET 

I TEMP. °F 

8 

63 

4944 

14.9 

66 

86 

7 

68 

4946 

14.8 

66 

88 

6 

76 

4948 

14.8 

66 

98 

5 

84 

4952 

14.8 

66 

100 


* "Thrust meter readings with initial load = 2.11 Lbs. 

Use Fig. (20) for calculations. Calculated readings 
for data in Table (31 } are as under. ' * 



SJffiLE 32 
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Calculated Data For Experiment of Determination Aid 
Of Performance Parameters Of A Gaseous Pocket Motor 
Using Hydrogen and Oxygen, At Different 0/P Ratios 
Por m = 0.012 Lb/Sec. 


0/F 

RADIO 

| P c 1 “ 

yEgf/ Cm 2 | Kgf * 

-L f 

tt — r 

j sp f 

|Kgf Sec/l 
I. . Kgm l 

°F 

I c* 

5 M/Sec. 
1 

I u 

% eq 

I M/Sec. 

\ 

8 

5,28 

0.766 

141 

1.365 

1010 

1385 

7 

5.63 

0.776 

142.5 

1.30 

1075 

1399 

6 

6.175 

0.790 

145 

1.208 

1180 

1421 

5 

6.72 

0.808 

148.3 

1.136 

1282 

1455 



APPENDIX B 


DJIA, COMPUTER PROGRAMME AND CALCULATED 
READINGS POE THEORETICAL APPROACH 


B. 1 Determination Of Combustion Temperature 

and Composition of Product Gas 

. B.2 Pro sen Plow Approach 

B.3 Equilibrium Plow App loach 
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B.1 Determination Of Combustion Temperature 
and Composition Of Product Gas 


Numerical solution of the equation (2.1.27) Section II is 
to be obtained. , To solve the set, equations are so arranged that 
the diagonal elements of the resultant matri:: are non zero i.e. 




**3 

^4 

*5 

r. c 

^7 

T»« 

*b 

o 

C 

y\ t 

•V>2. 


0 

2* 5 

0 

O 

^8 

o 

C-bW 

Cbic 

O 

o 

0 

T>3 

o 

Q 

2 h e 

0 

o 


(■< » 
c«tlc 

Q 

-l 

1 

0 

Vz 

0 

O 

0 

o 

0 

0 

O 

1 

0 

0 

-t 

M z 

0 

o 

0 

0 

0 

© 

o 

o 

1 

o 

-% 

-Vz 

0 

0 

0 

O 

o 

o 

0 

0 

-% 

0 

O 

1 

0 

0 

o 

0 

p 

o 

0 

0 

\ 

0 

o 

1 

o 

o 

o 

t 

6 

0 

o 

'O 

-V 2 

o 

o 

1 

0 

o 


2n, Q 2*4 Q O n 7 o o ^" qAi 


Calc 




A 

411^5 
hi* *V 
A J>> 9*j 

AT 




CcA w^ c ^ 

|. (%Z2toS. ) 

i f it TtuC N 
™ ' Ri c«ilc / 

*»cfcSff) 

Kfc> 

rw) ca | C 


(b.i.i) 


is the final set. to be solved. 
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Solution of set (B.I.l) is obtained using the 
standard disc programme available at I.I.T. Ksipur, Computer 
Center Library. This solves the given set using Gauss - Jordon 
method. Following steps were carefully followed to achieve 
convergence for the solution. 

(1) Sum of the initial guess of n^ should be preferably" 
less than p Q (this is preferable). 

(2) The estimated temperature is held constant for 
several iterations, thereby allowing adjustment of 
partial pressures to some more reasonable values. 

(3) Variation of temperature then by steps from start- 
ing value of 5000 °K till some reasonable error 
limit is established between initial guess of 
temperature and' calculated value of temperature. 

(4) Oscillations in solution can be avoided by modifying 
the error limits mentioned in (3). 

Equilibrium temperature and canposition is to be de- 
termined at pressures of 50, 100, 150 and 200 Psia in combustion 
chamber for OER varying from 2, 6 and 8. £L1 necessary ther- 

modynamic data for all species are tahen from standard tables 
given in References (19) and (29). A suitable, temperature in- 
terval of 200 °K is chosen. Attached herewith is the developed 
Computer programme for system IBM 7044/1401, which will provide 
the method of solution. 



o o o o o o 


c IB JOB MAP 

CIBFTC MAIN 

COMBUSTION TEMP, A I 
AND COMBUST I UN PR 
RF AD I MG QF VAR I US 
STANDARD T A B L - 3 AV 
UNITS OF ENTHALPY 
UNITS OF SP H ; ATS 
DIMENSION AK iU**), 
SHOT I { 24), Hfl T 2 ( 2-0 , 
-HOT 8 { 24) , HOT 9 ( 24 ) , 
-CPT6I24), C P T 7 ( 2 4 ) * 


D COMP 13IT1DN Cf-TcRMINATION FCR SIVcN 0/ F RATIO 
S A I JR • F UP 112-02 j AS ROCK FT M0TOR„ 

V -,LU AS MOTT ION ED IN DIMENSION STATEMENT FROM 
si L A B L 

A KCA L/MOL-4 
A R CAL/OK-MOLEc 

A K *' ( , AK ^ ( c 1 - ) , AK4 { 24 ) ,AKS(24) ,AK6(24) » 

H.iTMiO ,H0T4( 24) ,H0T5(24) » HQT6 ( 24 ) , H0T7 ( 24 ) » 
CPTI ( A) ,CPT2( 24) ,CPT3 ( 24 ) , CPT4 (24 ) ,CPT5( 24) , 

C P T 8 ( 24 ) , C P T9 { 2*0 , T ( 24 ) , OFR (4 ) , A ( 11 , 12 ) , B ( 11 , 1 ) 


U(U)tP(Ii!»0(9) 

REAL N19) ,N J(9).NT 
DATA BTRUF/ 2« / 

READ 2, L 

2 FORMAT (12) 

READ 3,M 

3 FORMAT (II) 

READ 4, (TCI >, 1=1,24) 

4 FORMAT (BFIi.ai) 

READ 6, (AKI ( I ) , I = 1» 24) 

6 FORMAT ( 6(EI( «4, 2X) ) 

READ 6, ( AK2 ( I ) , 1=1,24) 
READ 8,(4K3(I),I=i,24) 

8 FORMAT (6(EH«5t IX) ) 

READ 8, (AK4(I ), 1=1*24) 
READ 8, (AK5( I ) ,1=1,24) 
READ 10, ( AK6( J ) » I =1 » 24 > 
IS, FORMAT ' CtlE14*8, IX)) 

>■ READ 12, ( HOT If I ) ,1 = 1,24) 
12 FORMAT (8F1C« 4) 

READ 12,( H0T2( I ), 1=1,24) 
READ 12, ( H3T3U ) ,1=1,24) 
READ 12, ( HQT4( I ) ,1=1,24) 
READ 1 2 , ( HOT 5 ( I), 1=1, 24) 
READ 12* ( H0T6( I), 1=1, 24) 
READ 12* ( HOTTI I ) ,1=1,24) 
READ- 12, ( H0T8( I ) , 1 = 1 ,24 ) 
READ 12, (H0T9(I), 1=1*24) 
READ 14, { CPTI (I), 1=1, 24) 
READ 14,(CPT2< I) ,1=1,24) 
READ 14* f CPT3( I ) ,1=1,24) 
READ U»{GPT4m, 1 = 1,24) 
14 FORMAT (8Fie*3) 

16 FORMAT (8FU?*4) 

READ 14,(CPT5C1), 1=1,24) 
READ 16, ( CPT6I I ) # 1=1, 24 ) 
READ 14* C CPT?( I) *1=1,24) 
REA0l6i(GPT8( 10 , 1=1 ,24) 
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READ 16, (CPT9( I ) ,1=1,24) 

READ 18, ( OFM JJ ) , JJ=: ,4 ) 

18 FORMAT ( 4F1* „i) 

DO illi JJ = ., t n 
CALL FLUN (1 '■> ) 

INDCX= :■ 

K= . 

PREPARATIONS FOR MAKING CURRENT VARIABLES*, 

DO 853 1=1,9 
853 N ( I ) = o 25 
ANhQW = l«t,}0 
OPTRA =3*. A-l* 

STARTING TEMPARATURL IS ASSUMED TO BE 5 OKELVIN* 

AT =T ( 24) 

AAK1=AK1( 24) 

AAK2=AK2( 24) 

AAK3=AK3( 24) 

AAK4=AK4( 24) 

AAK5=AK5( 24) 

AAK6=AK6( 24) 

AHOTl=HOTl{ 24 ) 

AHQT2=H0T2(24) 

AH0T3 = H0T3< 24) 

AH0T4=H0T 4{ 24 ) 

AHOT5=HQT5f 24) 

AH0T6=H0T6( 24) 

AH077=H0T7t24) 

AH07 8 =HOT 8(24) 

■ ' A«8T9=H0T9* 24 ) ' 

A0m=CPTll24> 

ACPT2=CPT Zi 24) 

ACPT3=CPT3( 24 ) 

ACPT4=CPT4(24) 

AC PT5 = CPT 5( 24) 

ACPT6=CPT6( 24) 

ACPT7=CPT7( 24 ) 

ACPT8=CPT8( 24) 

AC PT9 = CPT 9( 24) 

FORMULATION OF CURRENT VARIABLES ENDS HERE* 

CALCULATION FOR VARIOUS ELEMENTS OF MATRIX ACI, J) IS STARTED* 
68 SU MN=N ( 1 ) +N ( 2 ) +N ( 3 ) +H (4 ) +N{5 ) +N( 6) +M( 7 ) +N< 8 )+N{ 9) 

AR1TR=AAK1 

AR2TR=AAK2 

AR3TRpAAK3 

AR4TR=AAK4 

AR5TR=AAK5 

AR6TR=AAK6 

ATRUE»U6*5i: /OFRl JJ ) )*2^ 

FAC0lM=(tl6*U0/0FR( J J ) > €v_ 

CTRU£“FACOI M*2» 0 
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AACALC = 2.*M( _ ) +N( 2 ) + Z* *N { *t ) + M { 7 ) 

BACALC=N( 1)+N(2)+N(?. H-2 *M( 5 ) + N { £) 

CACALC=N( 3) +2« *li ( U ) + '){ 9 ) 

ACALC=AACALC/ANECW 
BCALC = BAC AL C/AN ;GW 
CCALC=CAC ALC/ANEGW 
PTRA=CPTRA 
PCALC =SUMN 

AR1CAL = (NHI) >/(N(*)»(N(5) =5)) 

AR2CAL= t ( N ( 2 ) ) *{ N ( 4 ) •**(„ 5 ) ) / ( N ( 1 ) ) 

AR3CAL=<N(3) ) / < { N ( 6 ) -••*> 5 ) * ( N ( 5 ) ** « 5) ) 

AR4CAL=(M(7) )/<N(4)** «, .) 

AR5CAL=( Nt 8 ) )/(N(5) **■ <,;■) 

AR6CAL={M(9) }/(N(4» ) 

HACALC=AHOTI* i ( 1 )+AHDT2-*N (2 ) +AH0T3*N(3) + AH0T4*N (4 ) +AHQT5»N( 5 ) + 
1AH0T6*N ( 6 ) + AHQT7*N ( 7 ) +AHQT 8*N( 6 ) +AHQT9*N ( 9) 

HCALC=HACALC/AN r :CW 

HTRUE= ( IS*": /OFR { J J ) )«• { AH0T4 ) +4HCT i+F ACOI M*AHQT6 
SCPNT = ACPT1*N( 1 )+ACPT2*N( 2 )+ACPT3*N(3) +ACPT4*Nt 4)+ACPT5*N( 5 )+ 
,.ACPT6*N { 6 )+ACPT7*,M{ 7 ) + ACP T8* N( 3 ) +ACPT9*N ( 9 ) 

C PREREQUISITE CALCULATIONS FOR ELEMENTS OF MATRIX ENDS HERE® 

C ACTUAL CALCULATIONS OF THE ELEMENTS OF MATRIX STARTS. 

DO 150, LL=1 » li 
DO 15 , ■ M M= 1 , 12 
A(LL t MM) = i*.»A‘ 

1504 CONTINUE 

At 1,1 ) =N{ 1) 

At 1* 2 ) =N( 2 ) 

At 1 » 3 ) = N( 3) 

At 1*4) =N( 4) 

At 1*5) =N( 5) 

At I » 6 ) =N{ 6 ) 

At 1,7)=N( 7) 

At 1*8) =Nt 8) 

At 1 » 9 ) =N ( 9 ) 

At 1,12 ) = ( PCALC )*ALOG(PTRA/PCALC) 

At 2* 1 ) = N( 1) 

At 2, 2 ) *=N( 2 ) 

At 2» 3 ) =Nf 3) 

At 2,5)=2. *N( 5 ) 

At 2* 8) =N( 8) 

At 2, 10 ) =- ( BACALC ) 

At 2*12) = ( BACALC)*ALOG( B TRUE/ BCALC) 

At 3»3)=N( 3) 

At 3* 6 ) = 2. *N t 6 ) 

At 3* 9 ) =Nt 9) 

At3,l*}=MCACALC) 

At 3*12 > = ( CACALC ) *ALOG ( CTR UL/CCALC) 

At 4* 1 ) 

At 4, 2) -1.00 



A t 4, 4 ) =+» ! 'o 5 

A( 4, 12) =AL3G< AR27R/AR2CAL ) 

At 5,1) si. 

At 5,4 ) = ■" 1 « ■ ' * 

At 5,5)=- 0.5 

At 5,12)=AL3G( AR1TR/AR1C .L) 

At 6 » 3 ) =1. ... 

At 6, 5 )=“Do5 f . 

A(6,6)=~io5' 

At 6 , 1 2 ) = A L 3 G t A R 3 T R / 4 R 3 C , \ L ) 

A ( 7 » 4 ) =**0 « 5 
At 7, 7) =1.00 

At 7, 12) =ALOG< AR4TR/4R4CAL ) 

AC8,5)=-0.5G 

At 8,8 ) =!• U‘U 

At 8,12)=AL0G( AR5 7R/AR3CAL ) 

At 9,6)=* 0.5 
At 9,9) =1. 0 

At 9,12)=AL3G( 4R6TR/AP.6C AL ) 

At 1 » 1 ) =2 o * N ( j. ) 

At 1; , 2 ) =M { 2 ) 

At , 4 ) = 2 o # N ( 4 ) 

A t 1 . , 7 ) =N ( 7 ) 

At l.v , 10 ) = - ( AACALC ) 

At it » 12) = t AACALC )*AL0G ( ATRUfr / AC A 1C ) . 

IF ( INDEX. ECUl) GO TO 45 
DO 47 1=1, lit 
47 Bt 1,1 )=At I, 12) 

CALL MATINV (A,l( ,8, i, DET) 

DO 49 J=l,l 
49 Xt J) = B( J, 1) 

DO 880 1=1,9 

Pt I )=X( I ) +ALOG (N ( I ) ) 

N( I )=EXP( Pt I ) ) 

880 CONTINUE 

Pt 10)=Xt 10) +ALOG { ANriQ W ) 

AN£QW=EXP ( P ( 1 ) ) 

SNPPA=N tl )+N(2)+N{3)+N(4)+N(5)+N(6)+N(7)+N{8)+Nt9) 
RERRQR=t { SNPPA“'CPTR A) / ( CPTRA ) ) * 1 CU O 
IF ( ( ABSt PERROR ) )o Lie 5 o ! ) INDEX=\. 

GO TO 68 

45 At 11, 1 ) =AH3TI*N( 3 ) 

* A-t 11, 2 )=AH0T2*N ( 2 ) 

At II, 3 ) =A H3T3*N(3 ) 

At 11 , 4 ) =AH0T4*N ( 4 ) 

Ar:nt5>*AH0T5»N(5) 

At 11, 6 ) =AHOT 6*N 1 6 ) 

At 11 , 7 ) =A H0T7*N(7) 

At li,8)=AH3T8*N(8) 

At 11,9 ) = AH0T9*N t 9 ) 
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A( 11j 1*j > =•■( HACALC ) 

A ( 1 1 , 1 1 ) = ( 5 C P NT* A T ) / 1 > < ' v ' 

A ( 11,12 ) = (HACaLC )*ALOG( HTf UE/HC A LC ) 

C THIS FINISHES THE CALCULATIONS OF ; LAMENTS OF REQUIRED MATRIX* 

DO 11 I=1,U 
11 B( 1 , 1 ) = A( I, i.2 ) 

CALL M ATI NV (A»L»8»M, Dr: T ) 

DO 23 J=ltU 
23 X( J)=B( J,l) 

C HERE STARTS CHECKING OF XU) AND INTERPOLATION IF ANY* 

CALL MAXM I N ( X , L , X M AX , X M I N } 

C CALCULATIONS OF NEW VALUES TU BE USED, STARTS HERE. 

DO 890 1=1,9 
P{ I)=X( I >+ALOG(NU) ) 

N ( I ) = FXP< P<I) ) 

890 CONTINUE 

P( 1 ) =X ( 1 ) +ALQG ( AN .Qnl) 

ANEQW = cXP ( P ( 1' )) 

P( II ) = X( iU +ALUG ( vr> 

CALCAT = EXP( P ( 1 3 ) ) 

PERRAT= ( ( AT- CALCAT) /AT ) * 

IF (Pf-RRAT) i* i. , fit., ... * 1 
xSK’A IF ( ABS ( P ERRAT ) . LFio 25 ) GO TO 66 
IF (ABS(PERRAT).LE.i5.U- ) GO TO 9< 

GO TO 911 
9U. INQEX*0 

DO 9t i 1 = 1,9 
N(I)=1.83 
9o 1 CONTINUE 
ANEQW = 1*0CS 
AT = AT +10* OS 
GO TO 929 
91- INDLX=D 

DO 911 1 = 1,9 
N( I) =1.00 
911 CONTINUE 
ANEQW«1.0i* 

AT =AT + 100* D{ 

GO TO 920 

1301 IF ( ABS (PERRAT )*LE* 1* 25) GO TO 66 
IF CABS(PERRAT)eLE.2 : . ) GO TO 93. 

GO TO 940 

930 INDEXES 

DO 931 1=1,9 
N( I ) =1*00 

931 CONTINUE 
ANEQW~1®&1> 

AT =AT"lir» Vi 
GO TO 920 

94tv INDEX=0 
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DO 941 1 = 1,9 
N ( I ) = U<>Lf 
941 CONTINUE 
ANEQW = 1« !» 

AT =A T** 1 L D . . 

9 2* 1= ( ( AT”4u* o 1 <)/2 Q g * ) + . v < 

C INTERPOLATION OF TABLE VALUES STARTS HERE® 

AAK1=AK1 { I) + ( { A T~ T ( I ) ) / { T ( I + . ) - T ( I ) ) } * ( A K? ( I+D-AKi. ( I ) ) 
AAK2=AKZ( I ) + ( ( AT-T { I) )/(T(I + i)*T (I )) ) * { AK2 ( I +1 ) -AK2 ( I ) ) 
AAK3=AK3( I) +< ( AT- T( I ) ) / { T { I + ; ) ■ T < I ) ) }* ( AK3 (I +1 ) -AKi ( I ) ) 
AAK4=AK4{ I > -M < A T— T ( I ) ) / (T ( I + i )» T (I ) ) ) * { AK4 ( I +1 ) -AK4 ( I ) ) 
AAK5=AK5( I) + ( (AT«TU ) )/( T ( I + ;.)*■ T (I ) ) ) * { AK5 ( I +1 ) -AK5 { I ) ) 
AAK6=AK6( I > + ( ( AT-T ( I ) )/<T{I+i)-T (I ) ) ) * ( AK6 ( I +1 ) -AKo ( I ) ) 
ACPTl=CPT 1 ( I)+( (AT-T (I ) )/(T( I+i )-T< I)) )*(CPT1( 1+1 )»CPT1(I ) ) 
AC PT2=CPT2( I) + ( (AT-T(I) ) / ( T( 1+1 )-T ( I) ) ) * (CPT2 { I +1 )- CPT2 < I ) ) 
ACPT3= CPTStI ) + { (AT T (1 ) ) / ( T ( I +1 )• T ( I ) ) )*{ CPT3( 1+1 )~CPT3 ( I ) ) 
ACPT 4=CPT •+{ I )+( ( AT-T { I ) )/(T< I + 1)-T(I) ) )* ( CPT4 ( I +1 ) -CPT4 ( I) ) 
AC PT5 =CPT 5{ I )+{ (AT-T( I ) ) / ( T C I +2 ) -T { I ) ) ) * < C PT5 ( I +1 )-C PT5 (I) ) 
ACPTo=CPTb( I ) + ( ( AT-T( I ) ) /(T( I +i ) ~ T ( I) ) ) *(CPT6( I +1)«*CPT6< I ) ) 
ACPT7=CPT?( I ) + ( ( AT“T ( I ) 1 / { T ( I +1 ) -T { I ) ) ) * < CPT7 ( I +1 )~CPT7 ( I ) ) 
ACPT8 = CPT3( I )+( ( AT--T ( I ) ) / ( T { I +1 ) -T U ) ) ) * ( CPT8 ( I +1 ) ~CPT8 ( I ) ) 
ACPT 9 = CPT 9( I )+( C AT-T C I ) )/(T( I+i)-T(I) ) > * (CPT9( I +1 )“CPT9 ( I ) ) 
CPT18 = ( ACPT..+CPTK I ) ) /2 
CPT2B = (ACPT 2+CPT 2 ( I ) ) /2.J 
CP T3B= ( AC PT 3+CPT3 ( I ) )/2 
CPT4B = ( ACPT4-t-CPT4(I ) )/2,t 
CPT5B=(ACPT5+CPT5( I > )/Z, 

CPT6B=( ACPT6+CPT6CI ) )/2 a ! 

CPT7B=( ACPT7+CPT7U ) )/2*V 
CPT8B = t ACPT 8+CPT fi (I ) ) /2« ; 

CPT9B= ( AC PT 9+CPT9 ( I ) ) /2,i 

AHQTl=HOT Jl( I ) + ( ( CPT*B* ( AT'“ T( I ) ) ) /l OU.OO) 

AH0T2=H0T2( I ) + ( ( CPT2B* ( AT T( I) ) ) /i Ou.UlH 
; AH0T3=H0T3( I) + ( t CPT3S*UT- Tl I ) ) ) 

AHQT4=H0T4( I ) + ( < CPT4B,*C AT- T( I ) ) ) /l 
AHQT5*HOT5( I ) + ( < CPT5B* UT-T ( I ) ) ) /li.Cwufct;) 

I AH0T6=H0T6( I ) + ( ( CPT6B* { AT* T ( I) ) ) /i> 00. Oft) 

| AH OT 7 = HOT 7( I ) + ( { CPT7B* ( AT T(I)))/1 01*. UL ) 

? AHOT8-HOT 8( I ) + ( ( CPTBB* ( AT T(I)l) /l bumUi, ) 

AHOT9=HOT9( I > + ( (CPT9B*( AT- T( I ) ) ) /l OU.SjO) 

C INTERPOLATION OF TABLE VALUES ENDS HERE. 

K=K+1 

IF (K.GT. 3vr ) GO TO 3-' C 
GO TO 68 
; 66 CONTINUE 

1- NT = ( NU)+N(2)+N(3) + N(4)+N(5) +N(6 )+N (7 ) +N( 8) +N( 9 ) J/ANEQW 

DO 85 a 1=1,9 
NJ(I )=N(I )/ANEQW 
5 XJ( I )=NJ( I) /NT 

fe' 850 CONTINUE 
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PRINT 851 

851 FORMAT (//IX,* COMBUSTION CHAMBER PRESSURE GFR *) 

PRINT 8 52 » C PTRA i C FR ( JJ ) 

852 FORMAT { IX, 2U4o 5) 

PRINT 95 a 

95\ FORMAT (* NO OF MOLES OF UTROGF N PRESENT AS IMPURITY *> 

PRINT 951 , FACOIM 
951 FORMAT (IX, E12« A ) 

PRINT 93 

93 FORMAT (* ABSOLUTE MAXIMUM ERROR ABSOLUTE MINIMUM ERROR *) 
PRINT 29, XMAX,XM IN 
29 FORMAT (IX, 21: lie 4 ) 

PRINT 110 

11C FORMAT (* NEW VALUES OF M ( I ) S *) 

PRINT 120. (N( I ),I=1 ,9) 

120 FORMAT ( 1 X » 9 , 1 2 u 4 ) 

PRINT i3 j 

13' FORMAT (* Nl'W VALU ■ "IF • ,Q. FORMULA WEIGHT *) 

PRINT 142 , AN \QW 
14' FORMAT (IX, .2 2. A) 

PRINT 15: 

15 FORMAT (* Nt to VALUE OF COMBUSTION TEMP, IN DEGREE KELVIN *) 
PRINT 160 , CAL CAT 

16 FORMAT ( IX, L:i4o5 ) 

PR I NT St * J 

8. FORMAT (* INITIAL GUESS FOR COMBUSTION TEMPARATURE *) 

PRINT 829 , AT 
8 2U FORMAT <1X,E14o5) 

PRINT 600 

6„,. FORMAT (» PER, ERROR IN COMBUSTION TEMPARATURE CALCULATION *) 
PRINT 6 Ilf, PER RAT 
510 FORMAT (IX, El2o4 ) 

PRINT 86k> 

861 FORMAT (* MOLL FRACTIONS FOR VARIOUS SPECIES AT COMB TEMP *) 
PRINT 87U, ( XJ ( I) ,1=1,9) 

870 FORMAT (1X,9E12o4) 

GO TO 30 L? 2 

■ :3 not.:: PRINT 3‘131 

3001 FORMAT (* NO OF ITERATIONS EXCEED THE SPECIFIED LIMIT * ) 

39u2 CONTINUE 
1111 CONTINUE 
STOP 
END 

CIBFTC MAXMIN 

SUBROUTINE MAXMI N ( X , L , XMAX, XMI N ) 

DIMENSION X(U) 

XMAX=ABS( X( 1) ) 

XMIN=XMAX 
DO 31 J*1»U 

IF ( ABStXt J ) )« GT.XMAX) XMAX=ABS ( X ( J) ) 
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IF { ABS<X ( J ) JoLT.XMIf!) AM I N= ABS ( X ( J ) ) 

31 CONTINUE 
RETURN 
END 

CFETCH MAT INV CCS999 

SUBROUTINE M4TI NV ( C » N , B* M , D S T c R R ) 

DIMENSION C U.1 t i 2 ) , B ( 1 1 , , ) , I P I VO T { .. 1 ) » I NDEX (11*2) 
DOUBLE PRECISION A( II , 1 * » , AM AX, T , SWAP , PI VOT 
DO 5 1=1, N 
DO 5 J=1,N 
5 A( I, J ) = C( I, J) 

81- IF { AMAX-DA BS ( A ( J , K ) ) ) 8 5 , O ,1* 1 
95 AMAX = DABS(A( J,K) ) 

32y DETER M=DE TERM* PI VOT /DABS ( D I VOT) 

DO 78R 1=1, N 
DO 78i J = l, U 
78 C( I> J)=A( I, J ) 

CDON - 
C ENTRY 

C NOW FOLLOWS THE N.C.SSAxY DATA CARDSo 


MATUt; 
MAT. v 
MAT. s: 
MATO*: 
MATL’.l 
MATt V| 
MAT f I 
MAT Vyll 
MAT i;i| 
MAT on 
MATO! | 


i 





V\4- t ‘'l, •*? 1 
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Calculated datas are tabulates as under. 

TABLE 33 

Theoretically Calculated Data Dor Equilibrium Temperature 
and Composition Of Product Gas For Hydrogen, Oxygen, 
Nitrogen Propellant System Chemical Rocket Motor, For 

P c = 3.40 Ata 0/F Ratio = 2.00 

T Calculated = 3390.0 °K A Calculated = 0.3742 


H 2 ° 

0.6393 

0.1880 

OH 

0.08411 

0.02474 

NO 

0.0004784 

0.0001407 

H 2 

1.9730 

0.5804 

°2 

0.002277 

0.0006698 

N 2 

0.002900 

0.0008528 

H 

0.6773 

0.1992 

0 

0.01987 

0.005844 

N 

0.0004574 

0.0001345 


SPECIES 


NUMBER OF MODES J MOLE FRACTION 
n. X. 

-i A 
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TABLE 34 


Theoretically Calculated Data Dor Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen, 
Nitrogen Propellent System Chemical Pocket Motor, Por 

P c = 5.40 Ata 0/p Ratio = 4.00 

T Calculated = 3607.4 °E A Calculated = 0.6758 

c 


SPECIES 

/ FJMBER op MODES T~ 

1 n i s 

MOHS PRACTION 

nr 

H 2° 

0.T085 

0.2672 

OH 

0.2648 

0.07789 

NO 

0.001755 

0.0005161 

H P 

1.2488 

0.3672 


°2 

0.03011 

0.008855 

-2 

0.002113 

0.0006215 

H 

0.8272 

0.2433 

0 

0.1162 

0.03417 

N 

0.0007764 

0.0002284 


TABIE 35 

Theoretically Calculated Data Dor Equilibrium Temperature 
and Composition Of Product Gas Dor Hydrogen, Oxygen, 
Nitrogen Propellant System Chemical Pocket Motor, Eor 

p » 3.40 Ata 0/P Patio = 6.00 


“ c 

T OaloulfAed 
c 

* 3329.85 °K 

A Calculated = 0.9210 

SPECIES 

\ NUMBER OP MODES 

I \ 

■JliOIE PRACTION 

lx. 

JL j 

II o 0 

0.9886 

0.2907 , 

NO ' 

0.002826 

0.0008310 

h 2 

0.8660 

0.2547 

«2 

0-.09643 

0.02836 

» 2 

0.001550 

0.0004559 

H 

0.7865 

0.2313 

0 

0.2409 

0.07085 

N 

0.0008284 

0.0002436 •. , 

OH 

0.4168 

0.1226 ' j 
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TABLE 36 

Theoretically Calculated Data Dor Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen, 
Nitrogen Propellant System Chemical Socket Motor, Por 


P c = 3.40 

At a 

0/P Ratio = 8.00 

T Calculated 
c 

= 3674.5 °K 

A Calculated = 1.1350 

SPECIES 

I NUMBER OP IDLES 

$ n. 

{ i 

ImODE FRACTION ~~ 

I L 

H 2 0 

1.0350 

0.304 5 

OH 

0. 5068 

0.1491 

NO 

0.003576 

0.001052 

H 2 

0.6426 

0. 1890 

°2 

0.1921 

0.05651 

*2 

0.001246 

0.0003665 

H 

0.6775 

0.1993 

0 

0.3401 

0.1000 

N 

0.0007427 

0.0002184 


TABLE 37 


Theoretically Calculated Data Por Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen, 
Nitrogen Propellant System Chemical Socket Motor, Por 


P c = 6.80 .Ata 


O/P Ratio = 2.00 

1' Calculated = 
c 

3532.5 °K • 

A Calculated = 0.7484 

SPECIES 

I NUMBER OP MOLES 

\ V. 

I " MOLE FRACTION 

I L 

H 2 0 

1.2640 

0.1860 

OH 

0.1833 

0.02696 

NO 

0.001040 

0.0001529 

h 2 

3.9620 

0.5827 • 

°2 

0.004282 

0.0006257 

*2 

0.005664 

0,0008329 , 

H 

1.3380 

0,1968 

0 

0,03946 

0.005802 

N 

0.001105 

,0.0001625 



TABLE 58 

and Composition ^f^Pimduet ^ as^to" ^4 u ^ librium temperature 
m rog en Propellait S y st2 

a 0/F Ratio = 4.00- 

A Calculated = I.5500 


P 0 = 6.80 
T c Calcul ated = 3742.8 °K 


SPECIES ~} W/mR OF HOLES - ' i MOLE EEACTIOH 
x n. I v 


h 2 0 

CE 

HO 

H 2 

°2 

N 2 

H 

0 

N 


Z. 

JL 


1.7750 

0.5711 

0.003795 

2.5220 

0.05837 

0.003955 

1.6320 

0.2324 

0.001790 


0.2611 

0.08389 

0.0005581 

0.3708 

0.008584 

0.0005817 

0.2400 

0.03417 

0.0002632 


TABLE 39 

E C uliib ^'« Te-permture 

ntogen Propellant System Chemist £££“*>2?^ 

At a - 


P c = 6.80 


T c Calculated = 

3796.8 °Z 

A Calculated = 1.846 

SPECIES 

vIUMEER OE MOLES 

J n i 

~1 

| MOLE E^AOTIOH 

1 i 

O 

CM 

W 

1.9560 

0.2876 

OH 

0.8623 

0.1268 

HO 

0.006 132 

0.0009017 

h 2 

1.7760 

0.2612 

°2 

0.1957 

0.02977 7 

K 2 

0.002820 

0.00Q4146 ■ 

H 

1.5230 

0,2239 

0 

0.4776 

0.07022 

H 

0.001770 

C. 00026 02 
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2 .ABLE 40 

Theoretically Calculated Data for Equilibrium Temperature 
and Composition Of Product Gas For Hydrogen, Oxygen, 
Nitrogen Pro pell ant System Chemical Rocket Motor, For 

P c = 6.80 Ata O/'P Ratio « 8.00 

T c Calculated = 3802.4 °K A Calculated = 2.260 


SPECIES 

JJSUMBBH OP MODES 

1 h 

I MOLE FRACTION “ 
1 *j 

, H 0 0 

C. 

2.0030 

0.2946 

OH 

1.0550 

0.15520 

NO 

0.007608 

0.001119 

h 2 

1.3200 

0.1942 

°2 

0.3849 

0.05661 


0.002174 

0.0003197 

H 

1.3400 

0.1971 

0 

0.6853 

0.1008 

N 

0.001608 

0.0002364 


TABLE 41 

I 

Theoretically Calculated Data For Equilibrium Temperature 
and Composition of Product Gas For Hydrogen, Oxygen, 
Nitrogen Propellant System Chemical!. Rocket Motor, Por 

P„ = 10.20 Ata 0/P Ratio = 2.00 

T Calculated = 3606.4 °K A Calculated = 1.120 

c 


SPECIES 

feUMBBR OP MOLES 

I n i 

“I MOLE FRACTION 
\ X . 

i 2 

H 2 c 

1.8700 

0.1833 

OH 

0.2919 

0.02862 

'NO 

0.001750 

0.0001716 

H 2 

5.9400 

0.5824 

<4 

0.007035 

0.0006897 .. 

'2 ' 

0.008235 

0.0008123 

H 

2.0150 

0.1976 

0 

0.06352 

0.>:06227 

N 

0.001848 

0.C00 1811 
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TABLE 42 

Theoretically Calculated Data Per Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen 
Nitrogen Propellant System Chemical Rocket Motor, Por 

P c = 10.20 Ata 0/P Ratio = 4.00 

T c Calculated = 3835.7 °K A Calculated = 2.0240 


SPECIES 

| NUMBER OP MODE" 

l ”-i 

r mode REACTION 

1 X- 

i 2 

O 

CM 

w 

2. 6200 

0.2577 

OH 

0.8852 

0.08678 

NO 

0.006026 

0.0005908 

H 2 

3.8160 

0.3741 

°2 

0.08855 

0.008681 

n 2 

0.005676 

0.0005564 

H 

2.4170 

0.2370 

0 

0.3509 

0.03440 

N 

0.002860 

0.0002804 


TABDE 43 


Theoretically Calculated Data Por Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen, 
Nitrogen Propellant System Chemical Rocket Motor, Por 


P c = 10.20 Ata 


0/P Ratio = 6.00 

T Calculated = 
c 

3902.4 °K 

A Calculated = 2.7720 

SIECIED 

FnUMEERW 

_J n i 

MODES jt MODE PR ACT ION 

I 

h 2 o . 

2.9110 

0.2853 

OH 

1.3580 

0.1332 

NO 

0.009429 

0.0009243 

H 2 

2*6750 

0.2622 

°2 

0.2808 

0.02752 

N 2 

0.004025 

0.0003946 

H 

2.2560 

0,2212 

0 

0.7040 

0.06901 

N 

0,002850 

• 0.0002794 



TABLE 44 

Theoretically Calculated Lata Lor Equilibrium Temperature 
azid Composition Of Product Gas Por Hydrogen, Oxygen, 
Nitrogen Propellent System Chemical Socket Motor, Por 

P c = 10.20 Ata 0/P Ratio 8.00 

T q Calculated = 3909.4 °K A Calculated = 3.3960 


SPECIES 

\ NUMBER OP MOLES 

l *1. 

X 1 

I MOLE FRACTION ~~ 

I X. 

H 2 0 

2.9890 

0.2930 

OH 

1.6560 

0.1624 

NO 

0.01170 

0.001147 

H 2 

1.9860 

0.1947 

°2 

0. 5606 

0.05496 

K 2 

0.003062 

0.0003002 

H 

1.9780 

0.1939 

0 

1.0130 

0.09936 

N 

0.002550 

0.0002500 

TABLE 45 


Theoretically Calculated Lata Por Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen,. Oxygen, 
Nitrogen Propellant System Chemical Rocket Motor, Por 


P c = 13.60 

Ata 


0/P Ratio = 2.00 

T Calculated 
c 

= 3586.8 °K 


A Calculated = 1.4970 

SPECIES 

1 NUMBER OP MOLES \ MOLE FRACTION 

h 2 ° 

2.4910 


0.1832 

OH 

0.3988 


0.02932 

NO 

0.002393 


0.0001759 

H 2 ' 

7.9610 


0.5853 

°2 

0,009021 


0.0006633 

■ n 2 

0.01098 


0.0108071 

H 

2.6420 


0.1943 

0 

0.08252 


0 .006068 

N 

0.002593 


0.0001906 
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TABLE 46 


Theoretically Calculated Lata Lor •Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen, 
Hitrogen Propellant System Chemical Rocket Motor, Por 

P c = 13.60 

At a 

0/P Ratio = 4.00 

T Calculated = 
c 

= 3927.4 0 K 

A Calculated = 2.7030 

SPECIES 

InOMHER OP MOLE S 
! n i 

l MOLE FRACTION 

e 2 0 

3.5000 

0.2573 

OH 

1.2070 

0,8878 

no 

0.008259 

0.0006073 

H 2 

5.1220 

0.3766 


0.11550 

0.008495 

** 

0.007404 

0.0005444 

II 

3.1760 

0.2335 

0 

0.4601 

0.03383 

n 

0.003965 

0.0002916 


TABLE 47 


Theoretically Calculated Lata Por Equilibrium Temperature 
and Composition Of Product Gas Por Hydrogen, Oxygen, 
nitrogen Propellant System Chemical Rocket Motor, Por 

p c = 13.60 

At a 

O/P Ratio = 6.00 

I Calouated = 
c 

3964.2 °K 

A Calculated = 3.6910 

SPECIES 

| NUMHEH OF MDIES 
4 n. 

l MOLE PRACTIOn 

I X. 

f - 1 

h 2 0 

3.8280 

0.2814 

OH 

1.8200 

0.1338 

no 

0.01308 

0.0009618 

h 2 

3.6130 

0.2656 

°2 

0.3858 

0.02836 

u 2 

0.005064 

0.0003723 

H 

2.9830 

0.2193 

0 

0.9490 

0.06977 

n 

0.003857 

0.0 302836 
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• TABLE 48 

Theoretically Calculated Data Bor Equilibrium Temperature 
and Composition Of Product Gas Bor Hydrogen, Oxygen, 
Nitrogen Propellent System Chemical Rocket Motor, Por 

p = 13.60 Ata C/P Ratio = 8.00 

c 

T Calculated = 4017.5 °K A Calculated = 4.5450 

0 


SPECIES 

\ NUMBER OP MOLES 
\ n. 

l MOLE REACTION 

l X. 

f 2 


H 2 ° 

4.0070 

0.2947 


CH 

2.2080 

0.1623 


NO 

0.01611 

0.001185 


H 2 

2.6910 

0.1979 


C 2 

0.7622 

0.05604 


T,T 

2 

0.00383 

0.0002859 


H 

2.5750 

0.1893 


0 

1.3340 

0.09808 


N 

0.003380 

0.0002485 



B.2 Prozen Plow Approach 

Solution to the problem discussed in Part 2.3, Section II 
is attempted here. Calculated readings mentioned in Tables (33) 
to (48) are to be taken as input data for these calculations. Ne- 
cessary the mo dynamic tables for species are to be taken from 
Reference (l9). A suitable temperature interval of 200 °K is 
chosen. Attached herewith is the developed computer programme for 
systfem ILL 7044/1401 which will provide the method of solutipn. 
Error limits were modified at times to reach convergence. All 
calculations are carried out for small (X Pe = Atmospheric 
pressure. 



CUTE WATFOR, REWIND 

TC MAIN 

PERFORMANCE ANALYSIS OF GAS ROCKET MOTOR USING H2-02 
PROPELLANTS WITH FROZEN FLOW TECHNIQUE. 

READING OF VARIOUS VALUES MENTIONED IN DIMENSION STATEMENT ARE 
FROM STANDARD TABLES AVAILABLES . 

UNITS OF ENTHALPY ARE KCAL/MOLE 
UNITS OF SP HEATS ARE CAL/OK MOLE. 

UNITS OF ENTROPY ARE CAL/OK-MOLE 

DI MENS I ON HOT 1 (24 ) * HOT 2 (24) * HOT 3 ( 24 ) , H0T4 ( 24 ) ,HCT5 { 24 ) , HQT6 ( 24 ) , 
1HOT7 ( 24 ) » HQT8( 24 ) ,H0T9 ( 24 ) ,C PT1 ( 24 ) ,CPT2 ( 24 ) ,CPT3 ( 24) , CPT4( 24) , 
1CPT5(24),CPT6(24),CPT7(24) ,CPT8 (24)»CPT9{ 24) , SO TCI (24) , S0TC2(24) 
DIMENSION SOTC3I24) ,S0TC4(24),S07C5( 24 ) ,S0TC6 ( 24 ) , SGTC7 ( 24) » 
1SOTC8 ( 24) ,S0TC9(24) ,H0TE1(24) , HOTE2 ( 24) , H0TE3 ( 24) ,H0TE4<24) , 
1HQTE5 ( 24) , H0TE6 ( 24 ) , H0TE7 ( 24 ) , H07E8 ( 24 ) , H0TE9C 24) , CPTE1 ( 24 ) , 
1CPTE2 ( 24) »C PTE3 ( 24) , CPTE4 < 24 ) , CP TE 5 { 24 > , CPTE6 ( 24} , CPTE7( 24 ) 
DIMENSION CPTE3(24) , CPTE9(24) ,S0TEI( 24) ,S0TE2(24) ,S0TE3(24) , 
IS0TE4(24) ,SOTE5(24) , SOTE6524) ,SQTE7(24) , SOTE8 ( 24 ) , S0TE9 (24 ) ,T(24) 
1ET(24),TC(4)»XJ(4,9) ,0FR(4>» W(9) 

READ 4, (T(I ) » I =1 » 24 ) 
t FORMAT { 8FIU.1 ) 

READ 12 » ( HQT1 ( I ),I=1,24) 

2 FORMAT (8F10.4) 

READ 12, (HQT21I), 1=1,24) 

READ 12, ( H0T3 ( I ), 1=1,24) 

.READ 12, ( H0T4 ( I ), 1 = 1,24 ) 

READ 12, ( HDT5 ( I ) ,1=1,24) 

READ 12 , { HOT 6(1), 1 = 1, 24) 

READ 12, (H0T7(I), 1=1,24) 

READ 12, ( H0T8(I ) ,1=1,24) 

READ 12, ( H0T9( I), 1 = 1, 24) 

READ 14, { CPTIU) ,1=1,24) 

V FORMAT (8F10.3) 

READ 14, (CPT2(I), 1=1,24) 

READ 14,(CPT3(I),I=1,24) 

READ 16, (CPT4U), 1-1,24) 
i FORMAT (8F10.4) 

READ 14, (CPT5(I ), 1=1,24) 

READ 16,(CPT6(I ),I=1,24) 

READ 14, <CPT7(I ), 1=1,24) 

READ 16,LCPT8( I),I = 1, ; 24) 

READ 16,(CPT9(I ),1=1,24) 

READ 18,(OFR(JJ) , JJ=l,4) 

) FORMAT (4F10.1) 

READ 20, ( SOTCK I), 1=1, 24) 

3 FORMAT ( 8 FI 0.4 ) 

READ 20 , ( S0TC2 (I), 1=1, 24) 

READ 20,( S0TC3U 1,1 = 1,24) 

READ 20, ( S0TC4( 13,1=1,24) 

READ 20 , ( S0TC5 (13*1=1,24) 
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READ 20, ( S0TC6{ 13,1=1,24) 

READ 20, ( SDTC?( I), 1=1, 24) 

READ 20, ( SDTC8( I ), 1=1,24) 

READ 20, ( $0TC9( I), 1=1, 24) 

READ 4, ( E T( I ) , I = 1 , 24) 

READ 12, IHOTEKI), 1 = 1,24) 

READ 12,(H0TE2(I ), 1=1,24) 

READ 12, ( H0TE3 ( I ),I=1,24) 

RE AO 12 » ( H0TE4 { I ), 1=1,24) 

READ 12 , { HDTE5 ( I ), 1=1,24) 

READ 12, (H0TE6U 3,1=1,24) 

READ 12, (H0TE7(I ), 1=1,24) 

READ 12, {H0TE8U 3,1=1,24) 

READ 12, (H0TE9(I 3,1=1,24) 

READ 14, (CPTEKI), 1=1,24) 

READ 14, <CPTE2( I ), 1=1,24) 

READ 14, (CPTE3U ), 1 = 1,24) 

READ 16, ( CPTE4( I ), 1=1,24) 

READ 14, ( CPTE5 { I ), 1 = 1,24) 

READ 16 , ( CP TE6 ( I ), 1=1,24) 

READ 14, (CPTEKI ), 1 = 1, 24) 

READ 16, ( CP TE8 ( I ), 1 = 1,24) 

READ 16,(CPTE9{ I ),I = 1,24) 

READ 20, ( SOTEK I ), 1=1,24) 

READ 20 , ( SDTE2 ( I ),I=1,24) 

READ 20, ( S0TE3{ 13,1=1,24) 

READ 20, (SQTE41I 3,1=1,24) 

READ 20, (S0TE5( 13,1 = 1,24) 

READ 20, ( SQTE6 ( I 3 , I =1 ,24 ) 

READ 20, ( SQTE7U 3,1=1,24) 

READ 20,(S0TE8(I J, 1=1,243 
READ 20, (S0TE9U 3,1=1,24) 

READ 22,(TC(I 3, 1=1, 4) * 

22 FORMAT (4F10.4) 

READ 24, ( (XJ(I, II), 11=1, 9), 1=1, 4) 

24 FORMAT (3(E10. 4,2X3 3 
READ 26, (W( 13,1=1,9) 

26 FORMAT (9F6.1) 

THIS FINISHES THE TABLE READING 

CALCULATION OF EXIT TEMPRATURE STARTS 

ATC=AT 

UGC=1 .9857 

CPTRA=3.40 

RHS=UGC*ALOG(CPTRA) 

DO 1111 JJ=1,4 
AT £=3000* 00 ■ 1 
119 ATC=TC( JJ 3 

I=( (ATC-400.00)/200.00)+l-00 
(ATE-400.00 )/200. 00) +1.00 

AS0TC1 = S0TC1( I ) + ( ( ATC-T ( 1 3 ) / ( T ( I +1 )-T (1)33*1 S0TC1 C I+li-SOTCK 1 3 3 
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ASQTC2 = S0TC2( I ) + ( (ATC-T ( I ) ) / ( T ( I +1 )-T ( I ) ) } * ( S0TC2 ( 1 + 1 )-S0TC2 ( I ) J 
AS0TC3=SQTC3( I )+{ (ATC-T ( I ) ) / ( T ( I +1 I— T ( I ) } ) * ( S0TC3 ( 1+1 I-S0TC3 II)) 
AS0TC4=SQTC4( I ) + { ( ATC-T ( I ) ) / { T ( I +1)~T(I ) ) ) *( S0TC4( 1+1 >-S0TC4( I ) ) 
AS0TC5=S0TC5( I >+((ATC-T(I )) / ( T { I +1 )-T( I ) ) ) * ( SOTC5 ( I +1 ) -SOTC5 { I ) ) 
AS0TC6=S0TC6 ( I ) + { (ATC-T { I > ) / { T(I +1 )-T( I ) ) )*(S0TC6( 1+1 ) -S0TC6 ( I ) ) 
ASOTC7=SOTC7( I ) + ( (ATC-T (I ) ) / (T< Ml )-T ( I ) } )*($0TC7( 1+1 J-S0TC7 ( I ) ) 
ASOTC8=SOTC8( I ) + ( {ATC-T (I ) } / ( T (I + 1 )-T ( I } ) ) * ( SOT C8( I +1 ) -SOTC 8{ I > ) 
ASQTC9=SOTC9( I )+( (ATC-T (I ) )/ { T ( I +1 )-T { I ) ) ) *( SOTC9{ 1 + 1 )-S0TC9{ I ) ) 
ASOTE 1 = S0 TE 1 { J )+ ( (ATE-ET (J)) /(ET ( J+l )-ET ( J ) > )*( SOTE1 ( J+l J-S0TE2 ( J ) 
1 ) 

AS0TE2=S0TE2( J>+( (ATE-ET( J) )/(ET ( J + l 5-ET { J ) ) }*( SOTE2 ( J+l > -SQTE2 ( J ) 
1) 

AS0TE3=S0TE3( J ) + ( (ATE-ET(J) ) / (ET ( J + l )-ET( J ) ) ) * { SOTE3 { J + l ) -SOTE3 ( J ) 
II 

AS0TE4=S0TE4( J )+( ( ATE-ET ( J) ) / ( ET { J + l )-ET{ J } ) ) *{ S0TE4 ( J + l )-S0TE4( J ) 
il 

ASOTE5=SOTE5( J>+ ( (ATE-ET( J) )/(ET { J + l )-ET ( J I) I * { SOTE5 ( J + l I-SOTE5 { J ) 
1) 

ASOTE6=SOTE6( J ) + ( ( A TE- ET ( J )) / ( ET (J + 1)-£T{J) ) )*( S0TE6 ( J+l J-SQTE6 ( J ) 
1 ) 

ASOTE7 = SOTE7( J ) + ( (ATE-ETt J)) /(ET (J + D-ETC J) ) )*( S0TE7 (J + l )-SOTE7( J ) 

i) 

AS0TE8=S0TE8( J )+( (ATE-ET( J) J/(ET l J+l )-ET( J ) I )*( S0TE8 ( J+l I-SGTE8 ( J ) 
x) 

AS0TE9=S0TE9( J) + ( ( ATE-ET ( J)) / ( ET ( J+l )-ET ( J )J) *( SGTE9( J+l 1-SOTE9 ( J I 
1 ) 

TERMI=(XJ (JJ,1 ) ) * { ASQTC1-ASOTE1 ) 

TERM2=(XJ (JJ,2) i * (ASOTC2-ASOTE2 ) 

TERM3=(XJ (JJ,3) )*{AS0TC3-AS0TE3) 

TERM4* ( XJ ( J J * 4 ) ) *( AS0TC4-AS0TE4 ) 

TERM5=(XJ (JJ,5) I *( ASQTC5-ASOTE5 ) 

TERM6=(XJ(JJ»6) > * ( AS0TC6-AS0TE6 ) 

TERM7= ( XJ (JJ»7n*(A$QTC7-ASOTE7) 

TERM8 = (XJ ( JJ,8n*(ASOTC8-ASOTE8 ) 

TERM9*(XJ (JJ,9) )*(AS0TC9-AS0TE9I 

SUMTS“TERMi+TERM2+TERM3+TERM4+TERM5+T£RM6+TERM7+T£RM8+TERM9 
PERROR* ( ( SUMTS-RHSJ/RHS >*100.00 
PRINT 600 

600 FORMAT (* PERCENTAGE ERROR *> 

PRINT 601 »PERROR 

601 FORMAT (1X,E14.5I 

IF (ABS(PERRORUGE.25.00) GO TO 120 
IF {ABS{PERRQR>.l£.2.0G> GO TO 122 
ATE=ATE-10.00 
GO TO 119 

120 ATE=ATE-100 .00 
GO TO 119 
122 PRINT 132 

132 FORMAT (* RHS SUMTS *) 

PRINT 134, RHS, SUMTS 


Sv 
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134 FORMAT UX»2E14.5) 

C CALCULATIONS FOR EXIT TEMPRATURE END 

C CALCULATION OF EXHAUST VELOCITY STARS 

AVMWT=XJ{ JJ,2>*W(23+XJ{ J J , 3 ) *W ( 3 ) +X J ( JJ,4)*W{4)+ 

IX J( JJ,5)*WC5)+XJ(JJ,6)*W{ 6)+XJ( JJ,7)*W(7)+Xj( J J , 8 > *W ( 8 ) +X J C JJ,9)* 
AW( 9) 

AC PT1 =CPT 1 ( I ) + ( ( ATC-T{ I ) >/(T ( 1 + 1 3-T( I ) ) )*(CPT1( 1 + 1 )-CPTl { I ) ) 

AC PT2=CPT 2( I) + ( ( ATC-T ( I ) ) /( T ( 1 + 1 )-T( I ) ) )*(CPT2( 1 + 1 J-CPT21 I ) 3 
ACPT3=CPT3( I 3 + { ( ATC-T ( I ) ) / CT { 1 + 1 )-T ( I J ) )* (CPT3 ( 1 + 1 )-CPT3< I) ) 
ACPT4=CPT4( I ) + ( ( ATC-T ( I ) ) / { T ( 1+1 )-T ( I ) ) }*(CPT4{ 1 + 1 )-CPT4t I ) ) 

AC PT5=CPT5( I J + ( ( ATC-T ( I ) ) / 1 T < I + l)-T( I ) ) )*(CPT5{ 1 + 1 )-CPT5 ( I ) ) 
ACPT6=CPT6( I3 + ( { ATC-T ( I ) )/(T( 1+1 }-T{ I ) ) )*(CPT6t 1+1 J-CPT6 ( I ) ) 
ACPT7-CPT7I 1 3 + H ATC-T ( I ))/(T(I+l }-T(I)) )*(CPT7< 1 + 1 J-CPT7 1 1 ) ) 
ACPT8=CPT8C I) + ( { ATC-T { I M /< T ( 1 + 1 )-T ( I )))* { CPT8< 1 + 1 )-CPT8{ I) ) 
ACPT9=CPT9{ I > + { ( ATC-T ( I ) ) / { T ( 1 + 1 )-T ( I ) ) )*(CPT9( 1+1 )-CPT9( I ) ) 

CPT1B=( ACPT1 + CPT1 ( I ) )/2.0 
CPT2B = (ACPT 2+CPT2 ( I ) J/2.0 
CPT3B= ( AC PT3+CPT3 ( 1 ) )/2.0 
CPT4B = ( AC PT4+CPT4 ( I ) ) /2 .0 
CPT5B=( ACPT5+CPT5( 1 ) ) /2.0 
CPT6B=(ACPT6+CPT6(I ) )/2.0 
CPT76={ACPT7+CPT7(I ) J/2.0 
CPT8B = { ACPT8 + CPT8U ) J/2.0 
CPT9B = (ACPT9+CPT9(I ) )/2.0 

AHOT 1=HGT 1 ( I ) + ( (CPT1B*( ATC-T < I) ) 3/ 1000. 00) 

AH0T2 = H0T2( I ) + ( ( CPT2B* ( ATC-T (I) > 1/1000.00) 

AHQT3=H0T3( I ) + { ( CPT3B* { ATC-T < 1 13 3/1000.00) 

AH0T4=H0T 4( I ) + < ( CPT4B* ( ATC-T ( I) ) 3/1000.00) 

AH0T5=H0T5< I)+( { CPT 5B* ( ATC-T ( 1 3 ) 3/1000.00) 

AHOT 6=H0T 6( I ) + ( (CPT6B*{ ATC-T { I ) ) 3/1000.00) 

AH0T7=H0T7( I )+ ( (CPT7B* ( ATC-T ( I) ) 3/1000.00) 

AH0T8=H0T8< I )+( ( CPT8B*( ATC-T < I > ) 3/1000.00) 

AHQT9=H0T 9{ I ) + { ( CPT9B* { ATC-T £ I) ) 3/1000.00) 

ACPTE1=CPTE1( J)+( (ATE-ETt J33/(£T(J+13-ET( J3 ) 3 * ( CPTE1 { J+l 3-CPTE1 C J ) 
I) 

ACPTE2«CPTE2( J) + ( (ATE-ET< J) )/(ET U + l 3-ET { J) ) ) *{ CPTE2I J+l 3-CPTE2 ( J ) 
1) 

ACPTE3=CP TE3 ( J ) + ( (ATE-ETt J ) > / (ET (J+13-ETU )) ) *( CPTE3 { J+l 1-CPTE3I J 3 
1) 

ACPT E4=CPTE4( J ) + ({ ATE-ET { J ) 3 / (ET t J + l 3-ET ( J ) ) > *1 CPTE4 (J+1 ) -CPTE4U ) 

1 ) 

AC PTE5=CPTE5 ( J )+ ( (ATE-ET ( J ) ) / ( ET C J+l )~ET( J 3 ) 3 *{ CPTE5C J+13-CPTE5 < J) 

A) 

ACPTE6=CPTE6{J)+( (ATE-ET ( J ) ) / ( ET ( J+l )-ET ( J ) 3 )*( CPTE6I J+11-CPTE6UH 
1) 

ACPTE7=CPTE7( J)+l (ATE-ET { J))/(ET ( J+l )-ET{ J ) 3 ) *( CPTE7 { J+13-CPTE7 ( J) 



U 

ACPTE8*CPTE8 ( J 3 + ( (ATE-ET ( J ))/( ET (J+l 3-ET ( J ) ) 3 *( CPTEBl J+l 3— CPTE8C J) 
13 

AC PTE9=*CPTE9( J ) + { (ATE-ET ( J H / ( ET { J + l )-ET< J ) 3 3 * ( CPTE9( J+13-CPTE9 < JJ 
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i.) 

CPTE1B={ACPTE1+CPTEU J) 5/2.0 
CPTE2B=(ACPTE2+CPTE2( J) 5/2.0 
CPTE3B=(ACPTE3+CPTE3(J) 5/2.0 
CPTE4B= (ACPTE4+CPTE4{ J 5 5/2.0 
CPT£5B= ( ACPTE5+CPTE 5 ( J ) 5/2.0 
CPTE6B=(ACPT£6+CPTE6( J 5 5/2.0 
CPTE7B=(ACPTE7+CPTE7( J 5 5/2.0 
CPTE8B= ( ACPTE8+CPTE8 ( J ) 5/2.0 
CPTE9B={ACPTE9+CPTE9( J) 5/2.0 

AHOTE1=HOTE1( J 5+ ( (CPTE1B* (ATE-ET {J 5 55/1000.00) 

AH0TE2=H0TE2( J> + ( ( C PTE28* ( AT E-ET (J 5 5 5/1000.00 ) 

AH0TE3=H0TE3{ J)+ ( (CPTE3B* ( ATE-ET (J) ) 1/1000.00) 

AH0TE4=H0TE4{ J 5+ { (CPTE4B* ( ATE-ET (J 5 5 5 /1000.00 ) 

AH0TE5=H0TE5 ( J 5 + ( (CPTE5B* (ATE-ET ( J 5 5 5 /1000.00 5 
AH0TE6=H0TE6{ J 5 + ( ( CPTE68* ( AT E-ET ( J 5 5 5 /1000. 00 ) 

AH0TE7=H0TE 7 ( J ) + ( ( CPTE7B* ( ATE-ET ( J ) 5 5 /1000. 00 5 
AH0TE8=H0TE8( J)+( (CPTE8B*( ATE-ET (J 5 5 5/1000.00 ) 

AH0TE9=H0TE9( J 5 + ( ( CP TE9 B* ( AT E-ET (J ) 5 5/1000.00) 

TERMHi= (X J ( J J » 1 ) )*( AH0T1-AH0TE1 ) 

T£RMH2= (X J( JJ,2 5 )*( AHOT 2-AH0TE2 5 
T£RMH3=(XJ( JJ, 3 5 5 * ( AHOT 3- A HOT £3 5 
TERMH4= (X J( JJ»4 ) )*( AHOT 4-AH0TE4 ) 

TERMH5=(XJ( J J » 5 5 )*( AHOT 5- A HOT E 5 ) 

TERMH6={XJ( JJ t 6) }*( AH0T6-AH0TE6) 

TERMH7=(XJ( J J » 7 ) )*( AHQT7-AH0TE7) 

TERMH8=(X J{ JJ, 8) )*( AHOT 8-AH0TE8 5 
TERMH9= (X J{ JJ,9) )*( AH0T9-AH0TE9) 

SUMTH=TER MH1+TERMH2+TERMH3+T ERMH4+TERMH5+TERMH6+TERMH7+TERMH8+ 
1TERMH9 

PERFORMANCE PARAMETER CALCULATIONS FOR SMALL ALPHA AND PEXIT=PATM 
MKS UNITS ARE USED HERE. 

VEXIT={ (SQRT(SUMTH) ) * ( SORT ( 8370. 0) 5 i / ( SORT ( AVMWT) ) *SQRT( 10QQ.0) 

£VEXIT=VEXIT 

SP I MP=VEX IT/9.81 

GSTAR=(AL0G(CPTRA5)/<AL0G(CPTRA)-AL0G( ( ATC/ATE) 15 
CGSTAR=(2.0/(GSTAR+1. 5 } ** ( ( GSTAR+1 . ) / { 2.0* ( GSTAR-1 .Q 5 > ) 
CHVELO=(1.0/CGSTAR) *SQRT t ( UGC*ATC 5 / (AVMWT*GSTAR H*SQRT< 4185.0 J 
THCOFF«EVEX IT/CHVELO 
PRINT 127 

127 FORMAT (//10X,* COMBUSTION CHAMBER PRESSURE OFR *5 
PRINT 131»CPTRA»CFR< JJ) 

131 FORMAT {1X.2E14.55 
PRINT 124 

124 FORMAT (* COMBUSTION TEMPARATURE EXIT TEMPARATURE *5 
PRINT 126, TC( JJ5.ATE 
126 FORMAT (1X.2E14.55 
PRINT 136 

136 FORMAT <* AVARAGE MOLECULAR WEIGHT IN GMS/MOLE *5 
PRINT 138.AVMWT 



138 FORMAT (IX, £14.5} 

PRINT 128 

12b FORMAT (* VEXIT EVEXIT SPIMP GST AR CHVEIQ THCOFF *) 
PRINT 130, VEXIT, EVEXIT, SPIMP, GSTAR,CHVELO, THCOFF 
1 3'J FORMAT ( 1 X, 6E14 .6 ) 

1111 CONTINUE 
STOP 
END 

CENTRY 

C NOW FOLLOWS THE NECESSARY DATA CARDS. 
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Necessary data for calculations is taken from 
the Tables (33)? (34), (35) and (36). Calculated data is 
tabulated as under. 


TABLE 49 

Theoretically Calculated Data For Performance Parameters 
For Hydrogen, Oxygen, Nitrogen Propellant System Chemical 
Docket Motor At Different G/F Ratios With Frozen Flow 



Approach, 

For p = 
*c 

3.40 At a. 



O/F RATIO | 

2 | 

4 1 

6 | 

8 

T 

c 

°K 

3393.00 

3607.40 

3629.90 

3674.50 

T 

e 

°K 

2470.00 

2640.00 

2640.00 

26 50.00 

— * 
W 

Cms/Mole 

.5.3095 

7.9776 

10.1400 

12.0460 

Y e 

M/Sec. 

3575.29 

3006.66 

2699.93 

2533.54 

u 

eq 

M/Sec. 

3575.29 

3006 .66 

2699.93 

2533.54 

I 

sp 

Kgf Sec./ 
Kgm 

364.453 

306.490 . 

275.222 

258.261 

— * 
y 


1.35033 

1*34250 

1.35170 

1.36442 

c* 

M/Sec. 

3407.95 

2872.57 

2549.80 

2346.02 

C F 


1.04910 

1.04668 

1.05888 

1.07993 



Becessary data for calculations 


is taken from 


the Tables (37), (38), ( 39 ) and ( 40 ). Calculated data is 
tabulated as under. 


TABLE 


Theoretically Calculated Data Por Performance Parameters 

f 1 ) f ^f? enj Cxy ? en > Mtro S. en Propellent System Chemical 
Rocket victor At Different Cyp Ratios With Prozen Plow 
Approach, Por p = 6.80 At a. 


0/g MIC } 2 { 4 { 6 I 


T 

c 

°K 

L 

3532.50 

3742.80 

1 L 

3796.80 

3802.40 

T 

e 

°K 

2240.00 

2280.00 

2290.00 

2290.00 

— * 
W 

G-ms/Mole 

5.3037 

7.9672 

10.1650 

11.9970 

V 

G 

M/Sec. 

4226.13 

3686.07 

3322.35 

3069.96 

£ 

CD 

M/Sec. 

4226.13 

3686.07 

3322.35 

3069.96 

I 

sp 

K£f Sec./ 
Kgm 

430.798 

375.746 

338.670 

312.942 

— * 
Y 


1.31171 

1.34874 

1.35825 

1.35967 

c* 

M/Sec. 

3514.72 

2923.15 

2600.13 

2394.31 

c p 


1.20241 

1.26099 

1.27776 

1.28219 
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Hecessary data for calculations is taken from 
the Tables (4l),. (42), ( 43 ) and (44). Calculated data is 
tabulated as under. 

TABLE 51 

Theoretically Calculated Data Dor Performance P ar ameters 
Dor Hydrogen, Oxygen, Mtrogen Propellant System Chemical 
Rocket Motor At Different C'/P Ratios Vith Frozen Plow 



Approach, 

Por n = 

* c 

10.20 At a. 




0/P RATIO 1 

2 i 

f l 

4 

T 

1 

6 ~T 
6 f 

8 

T 

c 

°K 

3606,40 

3635.70 


3902.40 

3909.40 

T 

e 

°K 

2060.00 

2230.00 


2240.00 

2240.00 

— ■¥: 
w 

Gms/Mole 

5.3005 

7.9645 


10.1730' 

12.0130 

V 

e 

M/Sec. 

4611.94 

3863.10 


3490.50 

3325.95 

u 

eq 

M/Sec. 

4611.94 

3863.10 


3490.50 

3325.95 

I 

sp 

Kgf Sec./ 
Kgm 

470.127 

393.792 


. 355.810 

328.843 

— * 
y 


1.31775 

1.30468 


1.31411 

1.31544 

c* 

M/Soc. 

3546.64' 

2994.34 


2665.65 

2454.35 

n 

F 


1.30037 

1.29013 


1.30943 

1.31438 
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Necessary data for calculations is taken from 
the Tables (45), (46), ( 47 ) and (48). Calculated data is 
tabulated as under. 

TABLE 52 

Theoretically Calculated Data Por Performance Parameters 
Por Hydrogen, Oxygen, Nitrogen Propellent System Chemical 
Rocket Motor At Different O/p Ratios With Prozen Plow 


.Approach, 

Por p = 

0 

13.60 Ata. 



0/P RADIO | 

2 5 

4 “T 

4 l 

6 | 

8 

T ' °lv 

c 

3686.80 

3927.40 

3964.20 

4017.50 

T °K 

e 

2010.00 

2070.00 

2200.00 

2220.00 

W Cms/Mole 

5.3096 

7.9780 

10.1570 

12.0590 

M/Sec . 

4802.73 

4149.20 

3 5 96 . 30 

3349.23 

u M/Sec. 
eg. 

4802.73 

4149.20 

3596.30 

3349.23 


422.956 

366 . 596 

341.409 

— * 
y 

1.30279 

1.32515 

1.29133 

1,29409 

C* M/Sec.- 

3597.26 

3010.86 . 

2705.32 

2497.67 

0 

T? 

1.33511 

1.37808 

1.32934 

1.34094 


B.3 Equilibrium Blow Approach 


Solution of the problem discussed, in Part 2.4 , 
Section II is attempted here. Calculated data of equilibrium 
combustion temperature and various molefractions mentioned 
in Tables (33) to (48) are to be included here as data. Nece- 
ssary thermodynamic tables for species are to be taken from 
reference (l9). Suitable temperature interval of 200 °E is 
chosen. Special care is to be taken for assumption of initial 
guess of number of moles for various species. Over and above 
their sum, which should be less than c-xit pressure of one at- 
mosphere, they should be so chosen that MS constants of the 
set of equations should approximately come to zero. Error 
limits are suitably modified to reach convergence. Attached 
herewith is the computer programme developed for the computer 
system IBM 7044/1401. All calculations are carried out for 
small Cl and p. = Atmospheric Pressure. 
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CIBJCB MAP 

C IB FTC MAIN 

C PERFORMANCE ANALYSIS OF GAS ROCKET MOTOR USING H2-Q2 PROPELLANTS 

WITH EQUILIBRIUM FLOW TECHNIQUE. 

READINGS OF VARIOUS VALUES MENTIONED IN DIMENSION STATEMENTS ARE 
FROM THE STANDARD TABLES AVAILABLE. 

UNITS OF ENTHALPY ARE KCAL/MOLE. 

UNITS OF SP-HLAT ARE CAL/ OK MOLE. 

C UNITS OF ENTROPY 'ARE CAI./UK MOLE. 

DI MENS I ON HOT a ( 24 ) , HOT 2 ( 24 ) , HOT 3 (2 - ) , HGT4 ( 24 ) , HCT 5 f 24 ) ,H0T6(24) , 
iHoT7(24), HOTBf 24 ) ,HOT9 ( 24 ) , C PT* ( 24 ) ,CPT2(24) , CPT3 ( 24 ) , C PT4 ( 24 ) , 
aCPT 5 ( 24 ) » CPT6(24 ) »CPT7 ( 24} ,CPT6( 24 ! , CPT9{ ?4 5 , SO TCI ( 24) ,SuTC2(24) 
DIMENSION SOTC3I 24) , SOTC4 < 24 ) , SQTC5 ( 24 ) ,S0TC6(24 ) ,S0TC7 (24) , 
1SOTC8 ( 24) »SQTC9 ( 24 ) » HOT . i ( 24 ) , HOT E 2 ( 24 ) , HQTE3 ( 24 ) , H0TE4 ( 24 ) , 
1HOTC5 (24) ,HQT26( 24) * HOT cl ( 24 ) , HO Tr 8 ( 24 ) , NOTE 9 ( 2 A ) 

DIMENSION CPTEi (24) , CP Tl:2 ( 24 ) , CP TE J ( 24 ) ,CPTE4 ( 24 ) , CPTE5 ( 24 ) , 
iCPT66( 24) ,CPTt7( 24) ,CPTL8(24) ,CPTES(24) ,S0TE1(24) ,SOTE2 (24) , 
4.S0TE3 (24) »S0Tc4(24) , SOT: 5 (24 ) ,SGTci (24 ) ,SOT07( 24) , S0T£8(24) , 

LSOTf 9 (24) , AKEi ( 24) t HOF ( 9) , HOT RE ( 9 ) , T (24 > , £T ( 24 ) , TC ( 4 ) , QFR ( 4 ) 

DI Ms NS I ON W (9 ) ,X J(4, 9) , AO ( i ), 11 ) ,B;:( 10* 1) ,XE (10 ) ,PE( 10) * XJE (4,9 ) » 
iAVMWTCf 4) ,A VMWT2 (4 ) ,AMWR( 4 ) » AVMWT ( ) 

DI M.. 43 I ON A K- .2 1 :4 ) , AK£3 (24), AKc.4 (24 ) , AKE5 (24 ) ,AKE6(24) 

Rf.aL Nr (9 ) , NJ . ( '■) ) , M T L 

read •+ , ( t ( i i , i = i » _4 ) 

4 FORMAT ( BF _ . * > 

KsAD 6, ( A Kr: _ ( I ), 1 = 1,24) 

6 FORMAT (6 (c W.4, 2X) ) 

READ 6, (AK62U ) ,1=1,24) 

READ 8, (AKE3 ( I ) , 1 = 1,24) 

8 FORMAT (6 (Ell. 5, IX) ) 

READ 8, ( A KF4 ( I ) ,1 = 1,24) 

READ 8, ( A KE5 ( I ) , 1 = 1,24 ) 

READ 10, { AKE6U ) ,1 = 1,24) 

10 FORMAT (4(E!4.8,1X) ) 

READ 12, ( HOIK I ) ,1*1,24) 

12 FORMAT ( 8 F 1 0 . 4 ) 

READ 1 2 , { HO T 2 ( I ) ,1 = 1,24) 

READ 12, ( H0T3 ( I ) ,1=1,24) 

READ 12 , ( H0T4 ( I ) ,1=1,24) 

READ 12, ( HO T 5 1 I ) ,1 = 1,24) 

READ 1 2 , ( HO T 6 ( I ), 1 = 1,24) 

READ 12, ( HOT? { I ) ,1=1,24) 

READ 12, ( HOTS ( I ) ,1=1,24) 

READ 1 2 , ( HO T 9 { I), 1 = 1, 24) 

READ 14,(CPT1( I )»I=1»24) 

14 FORMAT (8F10.3) 

READ 14, ( CPT2 ( I ) , 1 = 1,24) 

READ 14, (CPT3( I )»I = lf 24) 

READ 16,(CPT4( I), 1=1,24) 

16 FORMAT ( 8 Fi 0.4 1 
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READ j.4,(CPT5(I),I = l,24) 
READ 16, { CP T 6 ( I ) ,1 = 1,24) 
READ 14, ( CPT7 ( I ) , 1 = 1,24) 
Kb AD J 6 , ( CP T8 ( 1 ) ,1=1,24) 
HEAD 16, ( C P T 9 { I ) ,I = _,24 ) 
READ i B » ( OF R ( J J ) , J J = , 4 ) 
18 FORMAT (4F1 >.„) 

READ 2 i , ( SO TC 1 ( I ) , 1 = 1 , 2 4 ) 
2U FORMAT (8 FI 0.4) 

READ 20, ( S0TC1 { I ),I=i,24> 
READ 20, ( S0TC3U )» 1=1,24) 
READ 20 , { SQTC4 ( I ) , I = 1 , 24) 
READ 20, ( S0TC5U ),I=1,24) 
READ 20, ( S0TC6U ),I = > , 2 h) 
READ 20, ( SDTC7 (I ) , I =1 , 24) 
READ 20, < S0TC6U ),I =1,24) 
READ 20, { S0TC9( I ) ,1=1,24) 


READ 

4, (£T( I ) 

,3 = 1,2'+) 

READ 

12, ( HOTE 

... ( I ) , I = 1 , 2 * ) 

READ 

1 R j ( HQ 1 l 

<!)»! = '. ,U ) 

R t A 0 

il,{ Hfil 

;■ ( x ) » 1 = r ) 

R c 1 •* 

8 , ( HO T .. 

i ( 1 ) f I ‘ :z i 1 u L \' ) 

R'i u Ij 

12 , ( HOT,. 

0(13,1 = :., 24 ) 

R c A 0 

L 2 , ( HO T L’ 

■A i ), l=i,24) 

RE AD 

1 2 » ( HO "1 . 

( 1 ) , I = a » 2 4 ) 

R l. A D 

i 2, ( HUT-. 

8(1 ) , I =1 , 24) 

REA'') 

a 2 , ( HO T E 

9 ( I ), I =a,24) 

READ 

±4, (CPTi; 

i( 1), I =1,24) 

READ 

14, ( CP It 

2(1 ) , 1=1,24) 

READ 

14, ( CP T E 

3 ( I ) , I =1 , 24 ) 

READ 

16, ( CP TE4 ( I 3,1=1,24) 


READ JL4, (CPTE5( 1 ) , I =1 , 24) 

READ 16, ( CPTE6 ( I }, 1=1,24) 

READ 14, ( CPTE7 { I 1,1 = 1,24) 

READ 16, ( CPTE8 { 1), 1=1,24) 

READ 16,(CPT£9( I), 1=1, 24) 

READ 20, « SOT Ei ( I 3, 1 = 1,24) 

READ 20, ( S0TE2 { I ), 1=1,24) 

READ 20, ( $OTO( I), 1=1,24) 

READ 20 , ( S0TE4 ( 13,1=1,24) 

READ 20,(S0TE6( I), 1=1, 24) 

READ 2u, ( SOT). 6 ( I 3,1=1,24) 

READ 2u , ( S0TE7 ( I ), 1=1,24) 

READ 20, ( SOT-.BU ), 1=1,24) 

READ 2 J, ( SO T 1 9 ( I >,I=i,24) 

READ 22,(TC(I) ,1 = 1,4) 

22 FORMAT (4F10-X ) 

READ 24, ( {XJ( 1,1 I), II =1,9), I =1,43 
24 FORMAT (3 (EU.4,2X) ) 


READ .26,1 HOF ( 1 ) , 1 = 1,9) 




0 0 0 0 0 0 
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4.6 FORMAT { 3Fi = j . 9 } 

READ 28 » ( HOTRt ( I ) , I = l , 9 ) 

28 FORMAT (9F8.4) 

READ 30, ( W( I ) ,1=1,9} 

3 FORMAT (9F6.1 > 

CP TR A = 6 .8 G 
tXIPR = j..OO 
UGC=1 • 985 7 

THIS FINISHES THE READING OF ALL NECESSARY DATA. 

CALCULATION FOR EXIT T E.MPARA TURC AND THEN FOR PERFORMANCE 
PARAMETERS STARTS HERE. 

TO START WITH WL WILL FIRST ASSUME THE EXIT TSMPARATURE AND THEN 
DETERMINE THE MOLE FRACTIONS OF VARIOUS SPECIES AT THAT EXIT 
TEMPARATURE AND EXIT PRLSSURc OF i ATA. 

DO 1111 JJ=i»4 
CALL FLUN (10000) 

AT £=2300. 00 

126 J= ( ( AT £-4 . . ) / 2 ' . ■')) + ,» . jO 

AAKE i = A KE j, ( J 5 + ( ( AT C - f.T ( J) ) / ( cT ( J +1 } - ET ( J ) ) ) * { AK El ( J + l >-AKE 1 ( J ) ) 
AAKc2= AKE 2 ( J 5 + ( ( AT-.-- T ( J ) ) / ( FT { J+i)-ET ( J 5 ) )*( AK E2 ( J +1 ) -AKE2 ( J) ) 
AAKFS = AKr 3( J 5 + ( (ATE- T ( J ) ) / ( c T ( J + 1 ) -t T < J ) > 5 * < AKE3 ( J + l J-AKE3 ( J 3 ) 
a A K ^ 4 = m K t **• ( J 5 + ( ( A T - - . T ( J ) )/(ET( j4i)-fcT( J) ) ) * ( AK E4 ( J + i ) -AK64 ( J 3 ) 
AAK L 5 = a K t 3 ( J ) + ( ( AT,-- IT ( J ) ) / ( cl { J + i ) -£ T ( J ) ) ) * ( AK E5 ( J + l )- AK£5 { J ) ) 
AA K = A KL^( J ) + ( ( ATE-LT (J) ) / ( i- T ( J + 1 ) -E T ( J ) ) ) * ( AK £6 ( J + 1 )- AKE6 ( J 3 ) 
KKK= J 

C INITIAL GUESS FCR VARIOUS VALUES OF NO OF MOLES. 

Nfc ( 1 ) =0.50 
NL: (2 5=0.1; 

IME ( 3 ) = J .0 .4 
NE (4 ) =J.l J 
NE (5) =0.10 

N6 (65=0.01 ' ' 

NE ( 7 ) = J .0 5 • 

NE ( 8 5 =0.05 
N£ (95=0.0001 
AN EQWE = 1 * 00 

128 SUMNE = Nfc( 1) +NE ( 23+NE ( 3)+Nt(4)+N£ (3 5+Nfc ( 6)+NE(7) +NE(83+NE(93 

ARE1TR=AAKEI 
AR E2TR = AA KE 2 
ARE3TR=AAK63 
: ARE4TR=AAKE4 

ARE5TR=AAKE5 
A ARE6TR=AAK£6 

ATRUSE = ( 16. OO/OFR { J J ) 3*2. iU 
BTRUEE=2.0D 

F ACOI M= ( (16. 00/OFR(JJ> 3+1.00 )*0.C01 
CT RUE 5= FA CO I M* 2. CO 

A A C A L E =2 . 00 * N E ( 1 > +N E ( 2 3 +2 . 09 *NE ( 4 3 + NE ( 7 3 
8ACALE=N£ ( i 3 +NE ( 2 ) + NiE ( 3 3 + 2 .00*N£ { 5 5 +N£ ( 8 3 
CACALE=NE (3 3+2.0C*M£(6)+N-(9) 
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AC ALc= AAC AL b / ANcC-Wi: 

BC ALt = BAC AL E/ANcGWE 
CCAL£=CACALt:/ANZCWH 
PTRAt=EXI PR 
PC ALE = SUMNE 

ARE 1CA= ( N E ( 1 ) ) / ( K t: ( 4 ) * ( 0 c ( 6 ) * * 0 . S 0 ) ) 

AREzCA= ( ( NE ( 2 3 )*(NE ( 4)**u.50 3 ) /( NE t i) ) 

A R fc 3 C A = ! N E ( 3 3 ) / ( ( N£ ( 6 ) ** > . 5 i ) * ( N E ( 5 ) ** . 5 • ) ) 
ARE4CA=(NE( 7) 3 / ( N E ( 4 3 # # ;. 5 j J 
ARE5CA= (NE{ 8 3 )/ C NE { 5 ) **0. 50) 

ARE6C A= (N E ( 9 3 ) / ( N t { 6 ) * *0. 5 •> ) 

C PREREQUISITE CALCUL AT I ON 3 END HERE. 

C CALCULATION FOR ELEMENTS :JF MATRIX STARTS HFRE. 

DO 70 11=1,10 
DO 70 KK= I, 11 
AE(II,KK3=0.Uv 
70 CONTINUE 

AE { 1 » 1 ) =N L ( 1 ) 

AE ( i , 2 ) =N ! r( t ) 

At ( 1, 3 )=Nt( i) 

At ( 1 , 4 ) =N ( ‘i ) 

A £ ( * , 3 ) =N : ! 0 ) 
kC ( 6 )=N '.( a ) 

Al ( i,7)=0t( 7 ) 

AE ( ) = lt( ) 

AE(i,9)=N£( 0) 

At (,. , j. i. ) = PC ALt* ALOG ( PTR AC /PCALfc ) 

A L ( A. i J. 5 = N c ( i, ) 

AE ( 2 , 2 ) =N E ( 2 ! 

A£ ( 2 , 3 ) =N E { 3 ) 

At (2, 5) =2 .0'J*Nfc ( 5 ) 

AE { 2 , 8 ) =NE( 8 ) 

AE ( 2, 10)=-( BACALE ) 

AE (2,11 3= (BACALE)*ALOG( BTRUEE/BCAL. 3 
Ac ( 3 » 3 )=NE( 3 3 
AE (3,63=2 «'Ju*NE ( 6 ) 

AE ( 3, 9)=NE( 9) 

AE (3 , 10 )= -( CACALE) 

AE (3, 11 3= (CACALE )*AL0G( CTRUEE/CC AL'c 3 
AE (4, 13 =-1.00 
AE (4 , 2 3 =1 .00 

A P / A . A J = + O . SO 

AE (4, 11 3 = Al OG ( ARE2TR/ ARE2CA ) 

AE ( 5,1 3=1 .00 
AE ( 5,4) =-1. 00 
AE ( 5 » 5 3 =-0. 50 

AE (5*113= ALOG (ARE ITR/AREiCA) 

AE (6,33=1-00 

AE (6, 5 3 =-0*50 

AE ( 6, 6 3 =-'Q. 50 ' • . 
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AE (6, U ) = ALOG( ARE'-iTR/AREjCA) 

At ( 7 , 4 ) =- 0. 50 
AE ( 7, 7 ) =1 .0 ) 

A£ ( 7 , 1 i ) = AL C Q ( A R E 4 T R / A R c 4 C A ) 

AE ( 8 » 5 ) =-Q. 50 
At ( 8 , 8 ) =1 . 9 „ 

At:' <8,11 ) = ALUG( ARE5TR/AR- 5CA ) 

AE (9,6)=-,. 50 
AE ( 9 , 9 ) =1 .00 

AE ( 9 , ]. i ) = AL 0G( ARS6TR / ARc 5CA ) 

AE (10,1 3=2. 00*NEU) 

AE ( 10 » 2 ) = N£ { 2 ) 

AE ( 10, 4) =2. 00*N£ (4) 

A E ( 1 0 » 7 ) = NE (7) 

AE ( 1 0 , 1 0 ) =- ( A A C A L c ) 

AE (10,11) =( AACALE )*ALOG ( ATRUtE/ ACAL E) 

DO 130 1=1, iO 
130 BE ( 1 , 1 ) =A £ ( I, iJ 

CALL MATI NV(A.;,_Oih „,COi ) 

DO j. 44 J=1,,.U 
144 :<E ( J) =B: ( J, i ) 

DO . 49 1 = j , 9 

PE ( I ) = X L ( 1 ) + /.LUG (Ot ( I ) ) 

NS { 1 ) =£XP (PL ( I ) ) 

149 CONTI MU. 

P!-(_. ) = y£' (l )+ALCG( A;JtQWc ) 

AN EQW 2= EX P ( Pc ( io ) ) 

C HE RE STARTS THE CHECKING FOR PP ADJUSTMENT. 

SNc PPA = NE { 1 )+Nfe ( 2 )+NE( 3 )+N£{ 4 )+NE{ 5 ) +NE (6 )+N£ (7 3+NE (8 1 +NE (91 
PERROE=( ( SNt PP A-EX I PR ) / EX 1 PR ) *10 C • 

IF ( ABS ( PERROE ) .LE.1.00) GO TO 162 
GO TO 128 

162 NT t = SN£PP A/ ANEQWE 
DO 166 1=1,9 
NJ E ( I )=NE (I ) / ANEC W£ 

XJE( JJ,I) =NJ£( I )/NTE 
166 CONTINUE 

AVMWTC ( JJ )=XJ ( JJ,1) *W(1 ) + XJ( JJ,2 )*W ( 2 )+X J < J J , 3) *W( 3 3+X J ( J J ,4 ) *W(4 ) 
i+X J ( J J » 5 ) *W { 5 ) + X J< J J,6>*W(6) +XJ( JJ,7)*wm+XJ( JJ,8 )*W(8)+XJ( JJ,9>* 
xW( 9) 

AVMWTE ( JJ ) = XJE ( JJ,i )*W( 1)+XJE( JJ,2 >*W(2)+XJ£( J J ,3 ) *W( 3 3 + XJg { J J ,4 )* 
1W(4)+XJE( JJ ,5)*W (5) +XJG (J J»6)*W( €)+XJ£( J J ,7) *W { 7 ) +X JE { JJ,S)*W( 81+ 
IX JE { J J » 9) *W ( 9 ) 

AVMWT( JJ) =( AVMWTC(JJ) + AVMWTE ( J J ) 3/2.00 
AMWR{ JJ )= AVMWTC ( JJ) /AVMWTE(JJ) 

C INTERPOLATION FOR ENTROPY STARTS HERE . 

AT C *T C ( J J ) 

I=( (ATC-400. 003/200.00)+!. 00 

ASOTCl = SOTCi( I ) + ( C ATC-T ( I ) )/(T(I+l )-T(I )) ) * ( S0TC1 ( 1 + 1 3-S0TC1 (1) 3 
ASOTC2=SOTC2 ( I ) + ( (ATC-TtI )) / ( T ( I 4 1 )-T ( I ) > )* ( SQTC2 ( I +1 >-S0TC2{ I) 3 



ASOTC3 = SOTC3 ( 1 ) + ( { ATC-T C I ) ) / { T { I *1 }-T { I ) ) ) * { S0TC3 ( I + 1 )-S0TC3( I ) ) 
AS0TC4=SGTC4( I ) + { { ATC-T U ) ) / { T ( I +1 )-T ( I ) ) ) * ( S0TC4{ 1 + 1 )-S0TC4( I ) ) 
AS0TC5 = SDTC 5 ( n+((ATC-T{I))/{T(I+l)-T(I)))*{S 0 TC 5 ( 1+1 >-SQTC 5 C 1 ) ) 
AS0TC6=S0TC6 { I ) + { (ATC-T ( I ) ) / (T { I +i)-T { ] ) ) ) *( SQTC6 ( 1 + 1 )-S0TC6 ( I ) ) 
AS0TC7=S0TC7 ( 1 ) + ( (ATC-T ( I )) / ( T ( I +1 )-T ( I } ) }*( S0TC7( 1+1 )-S0TC7( 1 ) } 
ASlJTC8=S0TC 8 ( I ) + ( (ATC-T ( I ) ) / ( T ( I +1 )-T ( 1 ) ) )*(S0TC3( 1 + 1 )-S0TC8( I } ) 

AS 0TC9 = S0 TC 9 ( I ) + ( ( ATC-T ( I ) ) / ( T ( I +1 )-T ( I } ) ) * ( SCT C9 ( I + 1 ) - S0TC9 < I ) ) 

J = KKK 

ASOTEi =S0 T'c 1 ( J) + ( (AT’-kT ( J) 5 / (IT (J + l ) — t T ( J ) ) ) *( S0TE1 ( J + l)-SOTfcl( J ) 
i ) 

A$OTE2»SOTE2( J ) + ( ( A Ti;- 1 T ( J ) )/(cT ( J + 1 ) - cT ( J ) ) )*( SQTc2(J + i)-SOTc2(J) 

x) 

ASOTE3 = SOTE3(J) + ( ( A T !-.-£■ T < J ))/( fc'T { J + 1 ) - 1 T ( J ) ) > * ( SOT c3(J + l)-SOTE3(J ) 
i) 

ASOTE4=SOTE4( J)+( ( A Tt-E T ( J ) ) / { £T ( J+1)-ET( J ) ) ) *( S0Tc4 { J + i )-SOT£4( J ) 

11" - 

ASOTE5 = SOTE5< J)+ ( ( ATL : - tT ( J )) / ( ET { J + l )-ET ( J 3 ) ) * ( SOTE5 ( J+ 1 )-SGTE 5 ( J ) 
i) 

AS0TE6=S0TE 6 ( J 3 + ( ( A T :-LT ( J ) } / ( tT ( J + 1J-ET( J) ) )*( S0TE6( J + l }-S0Tt6 ( J ) 
1) 

ASOTE7=SOTF7( J )+( ( ATf-ET ( J )) / ( ET ( J + 1 )-ET ( J ) ) ) * < SOTE7 < J + 1 J-SOTE7 ( J ) 

i ) 

AS OT E 8 = SO Tc b ( J ) + ( ( AT .- b T ( J ) ) / ( ET ( J + l > -ET { J ) ) )*( S0T£8( J + l )-SGTE8 ( J ) 
*) 

ASOTlr 9 = S0TC 9 { J ) + { (ATE-ET(J))/(ET ( J + 1 J-ET ( J ) ) ) *( S0T£9( J+l )-SQTE9( J ) 

i ) 

CHECK FOR irXIT T EM PARATURE STARTS HERE. 

RHSr>UGC*ALQG(CPTRA) 

TL-RSi = XJ( JJ,i)*( ASOTCl-UGC*ALOG(XJ( JJ»1) ) )-AMWR { J J )*X JE ( J J , U * 
i(ASOTEi-UGC*ALDG(XJE(JJjjt) )> 

TERS2=XJ( JJ,2 )*( AS0TC2-UGC*AL0G(XJ( JJ»2)> )-AMWR ( J J ) *X J£ ( J J * 21* 
x(ASGT£2-UGC*ALQG (XJt( JJ,2 >> } 

TERS3=XJ( JJ »3 )*< ASOTC3-UGC*ALOGlXJ< JJ»3) ) J-AMWR ( J J 1*XJ IE ( J J , 3 ) * 
i.(ASOTE3-UGC*ALOG(XJE( JJ,3) )) 

TERS4=X J ( J J ,4 )*( AS0TC4— UGC*ALOG(XJ( JJ»4) ) )-AMWR { J J ) *X J£ ( J J , 4) * 
i(AS0TE4-UGC*ALGG (XJE( JJ,4) ) ) 

. T£RS5=XJ( JJ ,5 )*( A$0TC5-UGC*AL0GCXJ{ JJ,5) ) J-AMWR ( J J >*X JE ( JJr 5 »* 

. iiASQTE5-JGC*ALOG{XJt{ JJ ,5) ) ) 

TERS6=XJ( J J , 6 ) * ( ASQTC6-UGC*ALQG ( X J ( J J* 61 ) l-AMWR £ J J ) *X JE { J J » 6) * 
1(AS0TE6-UGC*AL0G (XJE( JJ,6> ) ) 

T£RS7»X J ( JJ » 7 ) * ( A S0TC7-UGC*ALQG (XJ(JJ»7)> )-AMWR ( JJ )*XJE ( JJ,7)* 

1 { AS0TE7-U GC *ALOG ( X J fc ( J J » 7 ) ) ) 

TERS8 = XJ( JJ»S>*(AS0TC8-UGC*ALQG(XJf JJ,8) ) ) -AMWR t J J > *X JE { J J * 8) * 

I (ASOTE6-UGC*ALGG (X JE ( JJ,8 > ) 5 

TERS9=XJ( JJ * 9 ) * ( A SOTC9-UGC*ALGG( X J ( J J»9 J ) )-AMWR(JJ )*XJ€( JJ,9>* 

1( A$0TE9-UGC*ALGG (XJc ( J J *9} ) ) 

SUMS E = TER SI +T uRS2 + TERS3 + r _RS4+T£PS5+T£RS6+TERS7 +TERS8+TERS9 
PERR= ( ( SUMSE-RHSE )/RHSE )*10().0u 
PRINT 299 

299 FORMAT (* P ERR ATE *) 
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PRINT 300, PERK, ATP 
300 FORMAT (IX, 2El4.5 ) 

IF (ABS(PERR) .GF.60.00) GO TO 195 
IF (ABS(PERR) .L-.3. ./,•) GO TO L 75 
ATE=ATE + l0.iH> 

GO TO 126 

195 ATE=ATE+1 00.00 
GO TO 126 

175 PRINT 176 

176 FORMAT (//*0X.,* RHSr $UM$£ * ) 

PRINT 177 , R HS c , SUMS 1 

177 FORMAT (IX, 20*4.5) 

CALCULATION FUR EXIT TEMpARA TURE DETERMINATION ENDS HERE. 
CALCULATION FOR PERFORMANCE PARAMETERS STARTS HERE. 

INTERPOLATIONS NECESSARY ARE START ,D. 

ACPT1=CPT1< I)+( (ATC-T( I > )/(T( 1 + 1 )-T( I ) ) )*(CPT1( 1 + 1 3-CPTK I 5 3 
ACPT2=CPT2( I ) + ( ( ATC-T ( I 3 ) / ( T ( I +1 )-T ( I ) ) )*(CPT2( 1 + 1 )-CPT2 ( I ) ) 
ACPT3=CPT3( I 3 + < ( ATC-T ( I )}/(T(I + i 3-T ( I 3 ) 3 *( C PT3 ( 1 + 1 )-CPT3 ( I ) 3 
AC PT4=CPT4( I) + ( (ATC-T ( I) ) / ( T ( 1 + 1 )-T ( I ) ) )*{CPTA( 1 + 1 3-CPT41 I ) ) 
ACPT5=CPT5< 13 + ( ( ATC-T ( I 3 ) / ( T ( 1+1 )- T( I ) ) 3* (CPT 5 ( 1 + 1 1-CPT5 ( I ) ) 
ACPT6=CPT 6 ( I ) + ( ( ATC-T ( I > 3 / ( T (I + 1 3-T ( I 3 3 >*(CPT6{ 1 + 1 )-CPT6 ( I ) > 

ACPT7 = CPT7( I ) + ( ( ATC-T { I ) ) / ( T ( I + i )- I ( I 3 3 ) * ( CPT7 ( I + i ) -CPT7 ( I ) ) 
ACPT8=CPTY( I ) + ( ( ATC-T(I) 3 /(T( 1 + 1 )-T( I 3 3 )*(CPT3( I + i 3-CPTat I ) > 

AC PT9 = C PT 9 ( I ) + ( ( ATC-T ( I ) 3 / ( T ( I + 1 3-T ( I ) 3 ) * ( CPT 9 ( 1 + 1 )-CPT9 ( I 3 ) 

CP Ti 0= ( ACPTi+CPTKI ) 3/2.00 
CP TLB= ( AC PT 2 + CPT 2 ( I 3 3/2.00 
CPT3B = ( ACPT.; + CPT3(I 3 3/2-9 » 

CP T4B= { AC PT4+CPT4 ( I 3 3/2.0 J 
CPT5B= ( AC PT 5+CPT 5 ( I 3 3/2.00 
CPT6B=( ACPT6+CPT6(I 33/2.00 
CPT7B=(ACPT7+CPT7(I 3 3/2.0 1 
CP T8B= ( AC PT8 + CPT 8 ( I ) 3/2.00 

CPT9B=(ACPT9+CPT9(I 3 3 /2.00 - , , 

AH0T1 =HOT 1 ( I )+( (CPTxB*( ATC-T (13) 3/1000.00) 

AHOT2=HOT2( I )+( ( CPT 28* ( ATC-T ( I ) 3 3/ *000. 00 3 
AHOT3=HOT3( I ) + { { CPT3B* ( ATC-T ( 1 3 3 )/ 1(300. 00) 

AH0T4=H0T4( I ) + ( ( C PT 4B* { ATC-T ( I ) ) 3/ *000. 00 3 
AH0T5=H0T5( 1 ) + ( ( CPT58* ( ATC-T ( I 3 3 3 / 1000.003 
AHQT6-HQT 6( I 3 + ( ( CPT6B* { ATC-T { 1 3 ) 3/1000.00) 

AH0T7 = HQT 7( I ) + ( ( CPT 78* { ATC-T ( I ) ) 3/1000.00) 

AH0T8 = H0T 8 ( I ) + ( ( CPT8B* ( ATC-T (I ) ) 3/1000.00) 

AH0T9*H0T9< I ) + ( ( CPT98* ( ATC-T II ) > )/ 2000. 00 3 

ACPTEI*CPTE1( J)+UATE-ET{ Jli/tET ( J + l 3 -ETt 0 3 3 ) *( CPTEl C J+l I-CPTEl ( 0 > 

X } 

ACPTE2=CPTE2( J i+{{At£-ET{ J))/(ET ( J+l 3 -£T ( J 3 ) 3 *{ CPTE2 < J+l )-CPTE2( J > 

; 1 . i) ■ ■,.■■■ 

AC PTE 3® CP TE3 ( J ) + ( ( A T £ - E T ( J))/(ET ( J + l 3 -ET ( J 3 3 3 * < 0PTE3 ( J + l >-CPTE3( J ) 

: i) 1 ■ 

ACPTE4-CPTE41 J )+ { ( ATE-ET { J ) 3 / ( ET (J + l 3 -ETC J ) 3 3 *( CPTE4 { J+l )-CPTE4 ( J) 

1 ) ' ' . ' ' . ' 1 ‘ ' ' 



AC PTE 5=CP TE 5 { J )+ ( { AT -,-ET ( J ) ) / ( if T ( J + l ) - ii T ( J ) ) 5* ( CPT25( J + l i-C^TtS { J 5 

«» ) 

ACPTE6=CPTE6 ( J ) + ( ( AT. ,-£T( J 5 ) / ( c T ( J + i )-£T< J ) } ) * { CPTES ( J + i 5 -C PTE 6 ( J ) 

r \ 

4. f 

ACPTE7 = CP7: / { J 5+ ( ( AT C-cTf J 5 5 /{ rT ( J + l )--T( J ) ) )*( CPTE7( J + i 5-CPTE7 ( J) 

. 5 

ACPTbS=CPT6b ( J ) + ( { A'P - EI( J ) ) / ( cT ( J + i ) — iT ( J } J } *( CPT£3( J + i )-CPTS8< J ) 

i. ) 

AC PT 1 9=CP TF J( J 5 + ( ( AT --:i { J ) )/ ( C T (J + l )-uT( J ) ) }*( CPTc9( J + * 5-CPT29 { J 5 
*5 

CPT£1B=(ACPT£1+CPTE.( J ) )/?.. > 

CPTE2B=(ACPTfcL+CPT£2( J 5 5/2.90 
CPTt3B=(ACPTiri> + CP Th J{ J ) 5/2.00 
CP Tt4B= ( ACP Tc4 + CFT£ + ( J > 5/2.9+ 

CPTE5B=(ACPTE3+CPT£5< J 5 5/2.9 ; 

CPTE6B=(ACPT£6+CPT£b( J) 5/2.00 
CPTE7B= { ACPTcT+C P Tc M J ) 5/2.00 
CPTE8B=(ACPTfa + CPTEB( J) 5/2.0 + 

CPTE9B=(ACPTfc9+CPT£9( J ) 5/2.0 , 

AHGTE1 = HQTE1( J 5 + ( { CPTE1B* ( ATE-ET (J 5 5 5 /1000.00 ) 

AH0TE2 = H0Tt 2 ( J 5 + ( (C PT£2 B* ( ATC— tT { J 5 5 5 / 1 C00.C0 ) 

AH0TC3 = H0TE -A J5 + ( {CPTS3B* (ATE-tT (J) ) 5/iOOO.OO) 

AH0T£4=HQT£4{ J)+{ { C PTE4 B* < AT E-t T (J) ) 5 /i JO 0.00 J 
AH0TE5=H0TE5( J 5+ ( (CPT£5B*{ ATt-ET (J) ) 5/1000.00) 

AH GTE 6= HQ TE 6 ( J 5 + ( { C PTEfc 8* ( ATfc-tT (J 5 5 5 / 1000. DO) 

AHUTE7=H0T£ 7 ( J 5 + { ( C Pit 7 B* { AT S-tT ( J ) ) 5/1000.00 ) 

AH0Tt3=H0Tco( J 5 + ( ( C PTc8 8* { ATfc-fcT {05)5/1000.005 
AHQT£9 = H0T£9{ J 5 + { (CPT£9B*( ATt-ET (J) 5 J /1000.00) 

T£RHi = XJ{ JJ ,i )*{ EOF (I) + AHOT 1-HOT P£ (1 5 )-AHWR{ J J 5 *XJ£{ J J , 1 5 * { HOF Q > + 
i.AHQTbl-HQTRi.( i 5 5 

T £ RH 2 = X J ( J J j 2 5 * { HOF ( 2 5 + AHOT 2 -HOT RE { 2 5 )-AMWR< J J > *XJE { JJ,2) *( HOF { 25+ 
lAHOTE2-HOTRt{ 2 5 5 

T£RH3 = X J ( J J , 3)*(HOF(3)+AHOT3-HOTPE(3) 5-AMWR ( J J ) *X JE C J J , 3 5 * ( HOF 1 3 > + 
IAHQTE3-H0 TR fc ( 3 5 5 

T£RH4=XJ{ J J *4 >*{ FOF { 4 5 + AHQT4-B0T RE { 45 )-AMWR{ JJ) *XJE ( JJ ,4)*(H0F{4)+ 
iAH0TE4-H07R£{4) 5 

TERH5=XJ( JJ»8)*{F0F(5) + AHOT 5 -HOT RE ( 5 5 5-AMWRC J J) *X JE ( JJ ,5 5 * ( HOF { 5 ) + 
iAHOTgS-HQTRE (5) 5 

■TBIttt&kK J1 JJ ,6 5 *< HOF { 6 ) + AHOT 6 -HOT RE ( 6 5 >-AH WR { J J ) *X JE ( J J , 6) *{ HOF (65 + 
1 AHOTf fe-HO'TR ET 6 5 ) 

,f£RH7®XJC JJ *7 )* (HOF {7 J + AHOT7-HOT RE { 7 5 } -AMWR { J J) *X JE ( J J ,7 ) * { HOF ( 7}+ 
XAH0TE7- HOTR £ { 7 5 5 

TERH8=X J( JJ »8 5 * ( HOF { 8 5+AHGT8-HOT RE ( 3 ) 5-AMWR ( J J ) *XJE { J J , 8) * £ HOF { 8 ) + 
iAHOTES— HOTRE { B 5 ) 

TERH9=X J ( JJ » 9 } * { EOF ( 9 ) + AHOT 9— HOT RE { 9 ) )-AMWR ( J J) *X JE ( J J » 9 ) * { HOF { 9 ) + 

iAH0TE9-H0TRE ( 9 5 5 

SUMHE = TERH1+TERH2+T ERH3 +TERH4+TE RH5+TERH6+TERH7 +TERH8 + TERH9 
PERFORMANCE C ALCCLATI ON FOR SMALL ALPHA AND PEX IT=PATMA STARTS. 

MKS UNITS ARE USED HERE. 

V E X I T E - { ( SORT ( SUMHE 5 5 * { SORT ( 8370 .0 ) 5 / { SQRT ( AVMWT < J J ) 5 ) ) * ( SORT t 
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•Mb' A* 1-tf^ inJ I J ♦ VhJ J } 

LVtXIE = VEXIT£: 

SPIMPE=VEXlTL/9.8i 

GSTARc = ( A LOG ( CP TR a ) } / ( ALOG(C PTRA )-hLCG( ( ATC/ATt ) ) ) 

CG5TAE=(2 . ; /( GSTAkL : + i. ; ) )**( (GST AR_ + i . j )/( 2. J*( GSTARc-i . J ) ) ) 
CHV£L£= ( i . / COST A ' ) *SgRT ( ( UGC+ATC ) / ( AVMWT ( J J ) * 0 ST ARE ) )*SGRT ( Ax 85. G 

*) 

THCGFF = L-V -X J ‘ /CHV'L v 
PRINT 20 J 

20 j FORMAT (=5= C GMO t ST I fj- 1 CHAMBlR PRESSURE QFR *) 

PRINT 201, OPTRA, C FK ( JJ) 

201 FORMAT (IX, 2L14.5 ) 

PRINT 202 

2u2 FORMAT (* COMBUST ION T c PPARA I URc EXIT TE V PARATURE *) 

PRINT 203 , TC ( J J ) i AT w 
203 FORMAT { 1 X, 2G ^4 . 5 ) 

PRINT 400 

400 FORMAT (* MOLL FRAC1 1 OM OF VARiOLS SPECIES AT EXIT TEMPARATURE *> 
PRINT 401 , (X Jr! J J , I ) , I = _ ,9) 

401 FORMAT (IX, 9 .,9.5) 

PRINT 2 04 

2 ... 4 FORMAT {* AVMrtTL •* V M to 7 AMWR A VMWT *) 

PRINT 205 , A VMWT C { J J ) , A V ’Vil t: ( J J ) , AM wR ( J J } , AVMWT ( J J ) 

2,-5 FORMAT QX,4i. ,4.5) 

PRINT 2 6 

2 j 6 FORMAT (* Vt.XlT>- .VLXITb SPIMPE GST ARE CHVELE THCOFE *) 

PRINT 207, Vfc X i T _ , _V cX I £, SPIMPE, G STARE, CHVELE, THCOFE 
2,7 FORMAT ( IX, 6L..9 . 5 ) 

11 ii CONTINUE 
STOP 
END 

C FETCH M ATI NV CCS999 

SUBROUTINE MAT I iMV ( C , N, B , M , OETERM ) 

D1 PENSION C (10, 11), B( 10,1), I P I VO T ( 10 ) , I NDEX { 1 0, 2 ) 

DOUBLE PRECISION A ( 10 » 1 0 ) , AMAX, T ,SWAP , P 1V0T 
DO 5 1=1, N 
DO 5 J = i , N. 

5 A ( I, J »=C< I, J ) 

80 I F ( AM AX-DAB S ( A( J ,K ) ) ) 85 ,100, ICO 
95 AMAX=DABS (A ( J,K) ) 

32 J DE TER M=DETERM*PnOT/OA8S( PIVOT) 

00 780 1=1, N , 

DO 780 J=i, N 
780 C ( I , J ) =A( 1, J) 

COON6 
C ENTRY 

C NOW FOLLOWS THE NECESSARY DATA CARDS. 


MAT0002( 

MATC011C 

MATOOiif 

MAT0U13: 

MATS013* 

MAT0013< 

MAT0Q241 

MAT0027C 

MAT0O471 

MAT0077: 

MAT0077( 

MAT0U77S 



Necessary data for calculations is taken from 
the Tables (33), (34), (35) aid (36). Calculated data is 
tabulated as under. 


TABLE 33 

Theoretically Calculated Data. Tor Performance Parameters 

Por Hydrogen, Oxygen, Nitrogen Propellent System Chemical 

Rocket Motor At Different 0/P Ratios, With Equilibrium 

Plow Approach, p = 3. 40 Ata. 

c 


0/P RATIO | 

2 

1 1 

l 4 

\ 6 

l 6 

! — 

5 8 

T 

c 

C K 

3393.00 

3607.40 

3629.90 

3674.50 

T 

e 

°K 

3030.00 

3260.00 

3340.00 

3340.00 

— 

W 

Oms/iole 

5.4204 

8. '1573 

10.3730 

12.3230 

V 

e 

M/Sec 

3505.10 

2921.80 

2484.20 

2371.10 

u 

eq 

M/Sec 

3505.10 

2921.80 

2484.20 

2371.10 

I 

sp 

ICgf Secs/ 
Kgm 

357.30 

297.84 

253.24 

241.70 

— * 
y 


1.1019 

1.0902 

1.0730 

1.0846 

c* 

M/Sec. 

3627.40 

3060.90 

2738.90 

2518.20 


0.96626 


0.95455 


0.90704 


0.94160 
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Necessary data for calculations is taken from the 

tables (37), (38), ( 39 ) and (40). Calculated data is tabu- 
lated as under* 


TABLE 54 


Theoretically Calculated Data Por Performance Paramc-ters 
Por Hydrogen, Oxygen, nitrogen Propellant Slystem Ch emi cal 
Rocket Motor At Different 0/P Ratios, With Equilibrium 
Plow .Approach, p = S.tsO At a. 


0/P RATIO } 

L 

2 

T~T 

1 6 

l 6 

T 

8 

T c C K 

3532.50 

3742.80 

3796.80 


3802.40 

T e o K 

2950.00 

3220.00 

3220.00 


3310.00 

— * 

■W Grms/Mole 

5.4722 

0.2447 

10.5510 


12.4310 

v e M/Sec. 

4260.30 

3576.90 

3227.80 


2922.10 

U eq K ^ /Sec * 

4268.30 

3576.90 

3227.80 


2922. 10 

I Kgf Sec./ 
S P xr 

Kgm 

435.09 

364-v6l 

329.03 


297.87 

— * 
y 

1.1038 

1.0052 

1.0008 


1.0780 

C* ffi/Sec. 

3601,40 

3106,50 

276 9. 90 


2556 . 30 

°» 

1.1594 

1.1514 

1,1653 


1.1431 
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Necessaxy data for calculations is taken from the 
Tables (4l), (42), ( 43 ) end (44). Calculated data is tabu- 
lated as under. 


TABLE 55 

Theoretically Calculated Data Por Performance Parameters 

lor Hydrogen, Oxygen, Nitrogen Propellant System Chemical 

Rocket Motor At Different 0/P Ratios, With Equilibrium 

Plow ipproach, p = 1C. 20 Ata. 

0 


0/P RADIO | 



6 1 

0 

T 

c 

°K 

3606 #40 

3035.70 

3902.40 

3909.40 

T 

e 

°K 

2G00.C0 

3190.00 

3270.00 

3200.00 

— ** 






W 

Gms/Mole 

5.5050 

G.3C76 

10.6230 

12.5720 

v e 

M/Sec. 

4732.20 

3944.20 

3510.20 

3201.90 

u 

cq 

M/Sec. 

4732.20 

3944.20 

3510.20 

3201.90 

I 

sp 

Kgf Sec./ 
Kgm 

402.39 

402.06 

350.63 

334.55 

— * 
y 


1.1072 

1.0062 

1.0024 

1.0018 

c* 

M/Sec. 

3704.30 

3131.00 

2797.10 

2574.00 

°p 


1.2775 

1.2594 

1.2570 

1.2750 
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Necessary data for calculations is taken fron the 
Tables (45), (46), (47) and (40). Calculated data is tabu- 
lated as under. 

TABIE 56 

Theoretically Calculated Drta Tor Performance Parameters 
Per Hydrogen, Oxygen, Nitrogen Propellant System Chemical 
Socket Motor At Different 0/3? Patios, With Equilibrium 
Plow Approach, p = 15.60 At a 


0/P RATIO | 

2 r 

2 i 


6 r~ 

6 1 

0 

T c 

°K 

3686 .80 

3927.40 

3954.20 

4017.50 

T 

e 

°K 

2840.00 

3160.00 

3250.00 

3270 .00 

— * 
W 

&ms/Mole 

5.5291 

8.3655 

10.6850 

12.6400 

V 

e 

M/Sec. 

4979.20 

4203.70 

3758,80 

3419.40 

u 

eq 

M/Sec. 

4979.20 

4203.70 

3758.80 

3419.40 

I 

sp 

Kgf Sec./ 
Kgm 

507.56 

423.51 

383.16 

348.56 

— *- 
y 


1.1111 

1 .0909 

1.0824 

1.0856 

c* 

M/Sec. 

3732.50 

3153.10 

2811.10 

2599.00 

c p 

- 

1.3340 

1.3332 

1.3372 

1.3157 



